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Pilot  equivalent  time  delay  in  tin-  0 • Ks  loop;  decrees 
Total  open  loop  phase  ang 1 e of  t tie  a ► Fs  loop;  degrees 

n 

Phase  margin,  180”  + ) Y^,  (lwc);  degrees 
PSD  of  a,  (t) 

/p 

PSD  of  w^ft) 

PSD  of  A,.  ,(t ) 

t‘S 

PSD  of  0c(t) 

Imaginary  part  of  s;  rad/sec 

Break  frequency  of  the  turbulence  I’SD  model;  rad/sec 
Crossover  treqnency;  | Y ( J utc. ) | = I;  rad/sec 

Natural  frequency  ot  the  feel  system  modal  response;  rad/sei 

l'lie  minimum  frequency  for  which  l, dm)  = -1/N(m)  or  the 
minimum  frequency  for  which  i)>(]m)  = - 180”;  rad/sec 

Steady-state  frequency  of  fully-developed  P10;  rad/sec 

Resonant  frequency  of  pitch  attitude  closed  loop  pilot- 
vehicle  system;  rad/sec 

Undamped  natural  frequency  of  the  stick-1  ixed,  short -per i od 
mode;  rad/sec 

Undamped  natural  frequency  ot  the  stick-free,  short-period 
mode;  rad/sec 

I'.qu  iva  lent  natural  frequency  ot  the  N*'?,p(s)  polynomial; 
rad/sec  *' 


sec 


xv  l i 


T r ans  1 e r Kune  t i on  No  t a t i on 


(1/T) 
l t . u>  ] 


Short-hand  notation  for  [s  + (1/T)] 

Short-hand  notation  for  (a*  + 2 to  a 4 u>‘J) 

Numerator  polynomial  of  the  Xj (h)  transfer  function 


Mlsoel laneous  Notat Ion 


>1,0  ► Fs 


0 ► K„ 


The  closed  loop,  piloted  control  of  q anil  0 with  f' 
considered  as  the  pilot's  output 


The  closed  loop,  piloted  control  of  a_,  with  Fs  considered 
as  the  pilot's  output  *' 
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A.  BACKGROUND 

The  pilot- induced  oscillation  (PIO)  phenomenon  has  long  been  one  of  the 
most  fascinating  and  baffling  puzzles  in  manual  control  technology.  The  PIO 
problem  is  complex,  surrounded  by  controversy,  and  has  been  poorly  understood. 
It  is  a subject  that  is  very  important  to  understand  because  of  the  potential 
impact  a PIO  problem  can  have  on  an  aircraft's  development.  Historically,  PIO 
problems  have  tended  to  first  appear  in  the  final  stages  of  flight  test  and 
evaluation  and  are,  therefore,  very  expensive  to  correct.  There  has  been  no 
way  to  predict  PIO  encounter;  it  has  not  been  clear  whether  simulation  is  a 
viable  way  to  diagnose  PIO  tendencies. 

It  is  not  known  when  PIO  was  first  encountered  in  manned  flight.  The 
documentation  of  PIO  appears  to  have  begun  in  earnest  following  World  War  II 
at  a time  when  combat  aircraft  performance  was  being  extended  by  every  possible 
means  and  when  fully-powered  hydraulic  control  systems  were  being  tested  and 
introduced  into  service. 

A partial  list  of  United  States  aircraft  that  have  experienced  PIO- 
related  difficulties  is  given  in  Table  1.  In  some  cases,  the  PIO  resulted 
from  experimental  modifications  made  to  the  designated  aircraft  and,  therefore, 
is  not  necessarily  a characteristic  of  that  aircraft  as  it  is  normally  used. 
Many  other  undocumented  PIO  encounters  have  occurred.  It  is  probably  safe  to 
say  that  PIO  should  be  expected  to  occur  with  each  new  aircraft  until  the 
problem  becomes  sufficiently  understood  to  preclude  it  by  design. 

A bibliography  of  both  longitudinal  and  lateral-directional  PIO  is 
included  in  this  report.  The  reports  listed  are  representative  of  the  base  of 
data  and  information  available  on  the  subject.  No  attempt  was  made  to  list 
every  report  that  addresses  the  PIO  problem.  The  reports  are  listed  in  alpha- 
betical order  according  to  the  surname  of  the  first-listed  author.  A report 
listed  in  this  bibliography  will  be  referred  to  in  the  remainder  of  this 
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TABLE  1.  SUMMARY  OK  AIRCRAFT  KNOWN  TO  IIAVK 
ENCOUNTERED  PIO  (USA  Only) 


1. 

XP-42 

15. 

F- 1 06A 

•> 

P-63A-1 

lb. 

F-4C 

3. 

SB2C-1 

17. 

11-58 

4. 

F4U-4B 

18. 

1 

'-4 

> 

5. 

A 31) 

19. 

T-38 

6. 

C-97 

20. 

A4D-1/2 

7. 

F4H 

21. 

KC-135A 

8. 

A2F 

22. 

F-5A 

9. 

F-86D 

23. 

X- 1 5 

10. 

F-100C 

ro 

M2-F2 

11. 

F-101B 

25. 

YF-12 

12. 

F- 102A 

26. 

YF-1.6 

13. 

XF- 104 

27. 

YF-17* 

14. 

F-104B 

28. 

A- 7 

*As  simulated  with  the  NT-33A  only. 

section  bv  the  prefix  B(t‘or  example,  1112).  A separate  Reference  section  will 
list  principal  documents  supporting  the  analyses  of  this  report;  some  of  these 
will  also  he  found  in  the  Bibliography . 

The  1M0  problem  is  one  in  which  an  oscillation  of  the  aircraft  occurs 
(e.g.,  porpoising,  dutch-roll,  etc.)  that  is  difficult  or  impossible  for  the 
pilot  to  stop.  A central  characteristic  of  the  phenomenon  is  that  the  air- 
plane is  stable  both  stick-fixed  or  stick-free;  hence  the  name  p l lot  - induced 
osc l l lat ion. 

Opinion  varies  widely  on  the  causes  of  PIO.  The  dynamics  of  the  control 
system  have  been  implicated  by  numerous  investigations  as  a principal  cause 
for  riO  (e.g.,  114,  B5,  118,  1115,  1120,  B21,  1112,  11)4,  11)5,  B4b,  1157,  Bb5,  and 
1171).  Feel  or  control  system  nonlinearities  have  been  noted  as  important  to 
PIO  by  most  investigators;  among  these,  1115,  1154,  Bhh,  and  1171  are  ot  particu- 
lar interest.  Bobweights  have  been  identified  bv  several  investigators  as 
possible  PIO  initiators  (e.g.,  114,  111  5,  1)20,  1121,  11)8,  1)54.  and  1157).  Sur- 
prisingly few  citations  have  been  made  of  basic  airframe  dynamics  as  causes 
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for  PIO;  B46  is  one  of  the  few  references  to  suggest  that  airframe  properties 
are  at  fault.  In  HI'S,  the  role  of  the  airframe  is  considered  as  an  important 
part  of  the  overall  dynamic  system;  in  B22  a 1M0  criterion  is  proposed  that 
emphasizes  the  role  of  airframe  dynamics.  A few  studies  have  suggested  the 
importance  of  pilot  task  requirements  to  PTO  (B5,  B15,  B54,  Bbb,  and  B71). 
"Abrupt"  control  inputs  have  often  been  connected  with  the  onset  of  PIO 
(e.g.,  B15,  B27,  B54,  B66,  and  B71).  The  PIO  rating  scale  — shown  in 
Figure  3b  and  first  introduced  at  McDonnell  Aircraft  during  the  F-4  develop- 
ment --  emphasizes  the  significance  to  PIO  of  abrupt  control  (B21). 

The  piloting  cue  ttiat  is  most  central  to  PIO  is  a point  on  which  sharp 
disagreement  is  found  in  the  literature.  Normal  load  factor  response  has 
been  cited  direetlv  or  by  implication  by  numerous  investigators;  among  these, 
B4,  B“S,  B8,  and  B54  are  representative.  The  importance  of  load  factor  as  the 
primary  piloting  cue  is  discounted  in  B12  where  it  is  suggested  that  pitch 
attitude  is  the  primary  cue.  The  closed  loop  criterion  for  PIO  proposed  in 
B22  is  based  largely  on  pitch  attitude  tracking. 

One  notable  area  of  universal  agreement  (at  least  in  the  literature)  is 
the  importance  of  visualizing  the  pilot  as  a force  producing  dynamic  element. 

It  is  suggested  in  B15,  however,  that  the  pilot  should  also  be  considered  as 
a deflection  producing  servomechanism  in  analyzing  for  possible  PlO-prone 
a ircra f t -emit rol ler  comb inat ions . 

A lateral-directional  PIO,  known  to  exist  in  actual  flight,  has  been 
duplicated  in  a fixed-base  simulation  in  at  least  one  instance  (B47)  and 
possibly  others  as  well  (B73?).  However,  to  this  author’s  knowledge,  a 
longitudinal  PIO  has  never  been  created  in  a fixed-base  simulation,  although 
at  least  one  serious,  undocumented  attempt  to  do  so  was  made.  What  was 
called  a longitudinal  PIO  in  a fixed-base  simulation  was  obtained  by  Hirsch 
and  McCormick  (Hi1))  in  a fascinating  study:  this  was  a contrived  experiment, 
however,  in  which  rapid  and  non-physical,  otv-line  variations  of  1/Tp,  were 
made  in  order  t o catalyze  a sustained  system  oscillation — apparently  because 
a Pit'  could  not  be  initiated  otherwise!  The  non-existence  of  longitudinal  mode 
mode  I’ll'  In  tixed-base  simulations  is  regarded  bv  tills  author  as  a fact  of 
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some  Importance.  It  is  also  noteworthy  that  Fortenbaugh  (BJO)  failed  to 
create  PIO  in  a normal-acceleration-deprived  simulation  which  did  include  the 
pitch  rotational  acceleration  8. 


Ugl’l  . 


The  overall  state  of  the  art  of  PIO  understanding  may  be  best  inferred 
by  the  reader  after  he  has  studied  the  public  debate  by  Ashkenas  and  A'Harrah 
(BIJ,  B7,  and  B13).  This  valuable  record  suggests  that  there  are  many  ways 
in  which  the  nature  of  PIO  may  be  visualized  and  that  each  of  these  can 
contribute  to  our  ultimate  understanding  of  the  phenomenon. 

B.  THE  HANDLING  QUALITIES — Pit)  DICHOTOMY 

It  is  not  at  all  clear  whether  I’ll)  necessarily  results  from  basically 
poor  handling  qualities.  The  evidence  is  that  it  does  not.  That  is,  P10  has 
occurred  with  aircraft  that  have  passed  flight  test  certification  and  may 
therefore  be  presumed  to  have  generally  acceptable  handling  qualities  (the  T- 
58A  is  an  examp le- -Kef . 11. 

a i 
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On  the  other  hand,  if  an  aircraft  is  very  difficult  to  control  (Cooper- 
Harper  ratings  worse  than  about  S)  must  we  conclude  that  it  is  FlO-prone? 

Can  we  determine  t torn  flight  test  whether  such  an  aircraft  has  a 1’tO  problem 
or  whether  it  merely  has  poor  handling  qualities  in  the  control  of  attitude? 
Does  it  matter? 

l'here  are  two  complicating  factors  involved  in  attempting  to  answer  such 
questions.  First,  the  relation  between  the  pilot's  task  and  any  tendency  to 
excite  latent  F10  in  a flight  or  simulator  experiment  is  unknown  since  we 
know  neither  il  latent  PIO  exists  nor,  it  it  does,  how  it  may  be  excited. 
Second,  we  have  no  rule  for  discriminating  between  a control  problem  due  to 
poor  attitude  control  dynamics  and  one  due  to  PIO  tendencies  when  these  are 
encountered  in  t 1 ight  test  or  simulation.  Thus,  we  cannot  be  sure  that  we 

I are  testing  the  proper  control  situations  or  that  we  can  identify  a potential 

F10  coat  igur.it  ion  when  one  is  encountered.  This  is  a nontrivial  matter  since 
i it  mav  at  feet  the  character  ot  system  change  requirements  and  the  impact  that 

the  airciatt  might  have  it  delivered  to  the  training  or  service  common i t i es . 
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I'he  tendency  in  past  work  has  been  to  identify  any  pilot-aircraft 
oscillation  as  a 1’10  without  a great  ileal  of  discrimination.  Of  course,  the 
w ings-bend ing,  rivet-popping,  fully-developed  I’lO  is  easily  identified  when 
it  occurs  (by  definition).  This  is  not  the  usual  PIO  encounter,  however.  We 
need  a method  for  the  unambiguous  separation  ot  attitude  control  problems 
from  those  that  we  choose  to  call  t'10.  Clearly,  a workable  definition  for 
PIO  is  required. 

C.  REPORT  OUTLINE 


In  the  remainder  of  this  report,  only  the  longitudinal  mode,  short- 
period  Pit'  will  be  considered.  A formal  definition  for  PIO  will  be  proposed 
in  Section  11  which  will,  it  is  hoped,  permit  resolution  to  be  made  between 
the  PlO-prone  aircraft  and  the  one  that  merely  has  deficient  handling  qualities 
in  the  control  of  attitude.  A model  for  pilot  dynamics  will  be  presented  in 
Section  111  that  is  useful  for  the  analysis  of  PIO.  Connections  between  this 
model  and  the  more  familiar  ones  found  in  the  literature  of  pilot  dynamics 
will  be  d i will  be  dev.' ted  to  the  so-called 
synchronous  tracking  mode  of  pilot  behavior  that  has  been  linked  with  PIO  in 
past  works.  Section  IV  will  present  a physical  and  mathematical  theory  for 
PIO  that  is  consistent  with  available  data.  Two  classifications  of  PIO  will 
be  Introduced  in  an  attempt  to  account  for  the  mechanisms  by  which  a PIO  can 
be  catalyzed  in  actual  flight.  A compendium  of  rules  for  the  assessment  of 
PIO  will  be  summarized  in  Section  V as  an  aid  to  applications  of  the  PIO 
theory.  A number  of  numerical  examples  will  be  given  in  Section  VI  to  indi- 
cate that  known  PIO  problems  are  "predictable"  with  the  proposed  theory; 
equally  Important,  one  example  will  be  given  to  confirm  that  an  aircraft  not 
known  to  have  experienced  PIO  difficulties  is  indeed  excluded  as  a PIO  candi- 
date by  present  theory.  Section  VI 1 will  discuss  some  general  problem  areas 
of  importance  to  the  experimental  identification  of  PIO  from  the  viewpoint  of 
the  proposed  theory.  Finally,  Section  Vlll  will  summarize  what  has  been 
concluded  elsewhere  in  this  report. 


A derivation  ot  the  pilot-centered  normal  acceleration  transfer  function 
is  presented  in  the  Appendix.  Approximate  factors  for  the  acceleration 
transfer  function  are  also  derived  there. 
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SKCT10N  11 
A Dll' IN  ll'l  ON 

The  lack  ot  a suitable  del  inil  Lou  lor  I’lO  lias  been  a major  obstacle  it' 
nude  island  ill*.;  t hi1  phenomenon.  If  appears,  however,  that  a nontrivial  I’lO 
dot  mil  ion  can  only  be  tormulated  alter  the  physical  problem  is  understood 
with  a certain  degree  ol  sophist  icat ion.  We  arc  theretore  faced  with  the 
task  ot  developing  an  understanding  ol  a problem  we  cannot  vet  del ine. 

Needless  to  sav , this  author  went  down  many  blind  alloys  before  arriving  at 
the  thooret  ical  descript  ion  ot  I’ll'  that  t orris  the  main  body  ot  this  report  . 
based  i'u  the  results  ol  Sect  ion  IV,  then,  a consistent  I’lO  del  init  ion  can  be 
do v t sed  as  I o I lows  : 

A.  I’lO  DEFINITION 

I’lO  is  an  unwanted,  inadvertent  and  atypical  closed  loop  coupling 
between  a pilot  and  two  or  more  independent  response  variables  of  an 
a i re  rat  l . 

It  is  presumed  that  a fully-developed  I’ll'  will  either  tax  the  pilot’s 
ability  to  control  it--and  thereby  leave  little  or  no  margin  for  accomplish- 
ment ot  the  primary  piloting  t ask— or  it  will  b>  entirely  beyond  his  control 
capab i 1 t t les . In  the  latter  case,  the  only  success! ul  recovery  technique  is 
lor  the  pilot  to  remove  himself  from  the  control  loop,  either  by  clamping  or 
releasing  the  control. 

b.  COMMENTS 

This  definition  eliminates  from  the  I'll'  classification  all  those  single 
variable  feedback  control  systems  for  which  lightly  damped  closed  loop  oscilla- 
t Ions  are  intrinsic  (even  with  a high-gain  autopilot!)  but  which  are  not 
otherwise  susceptible  to  I'll'.  The  definition  of  I’ll)  as  "unwanted"  excludes  the 
"stick  pumping"  phenomenon  (Ref s. 2 and  ))  from  the  I’lO  category. 

The  "independent"  response  variables  referred  to  in  the  I’lO  definition  are 
intended  to  correspond  with  the  usual  aircraft  degrees  of  freedom,  their 
derivatives,  or  linear  combinations  ot  these.  A derivative  of  one  of  these 
variables  is  not , bv  definition,  independent  of  the  variable  itself. 
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By  this  definition  one  might  suspect  that  a certain  element  of  chance  is 
involved  in  PIO  encounter;  this  author  suspects  that  this  is,  indeed,  the 
case.  Thus,  we  must  test  for  PIO  with  great  care  to  he  certain  that  the 
conditions  necessary  to  sustain  it  are  satisfied. 

It  is  not  by  chance  that  this  definition  dictates  that  P10  requires  a 
level  of  piloting  task  that  is  more  complex  than  mere  control  of  attitude. 

It  is  contended  that  PIO  can  only  occur  for  those  circumstances  in  which  the 
pilot  is  very  concerned  about  precision  control  of  aircraft  path  or  path- 
related  acceleration  components. 

Past  occurrences  of  Pit)  have  been  with  aircraft  that  were  otherwise 
stable;  with  future  aircraft  employing  relaxed  static  margin,  PIO  could  occur 
following  control  system  failures  for  which  the  aircraft  would  be  unstable. 

Such  a control  problem  is  still  admitted  to  the  PIO  category  by  the  defini- 
tion here;  however,  for  such  aircraft  stability  could  not  be  restored  following 
PIO  encounter  by  clamping  or  releasing  the  controls. 
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SECTION  III 
A PILOT  MODEL 


A . BACKGROUND 


It  was  assumed  in  References  4 and  5 that  the  essential  piloting  cue  in 
PIO  is  pitch  attitude;  it  was  further  assumed  that  a sufficient  model  for 
pilot  dynamics  in  the  fully-developed  PIO  is  a pure  gain;  i.e.,  Y^(jw)  = 
such  that  the  stick  force  F (t)  is  just  K^O  (t).  No  pilot  dynamic  equali- 
zation (lead  or  lag)  or  time  delay  was  considered  necessary  for  analysis  of 
the  causes  of  a PIO  once  ij^  was  known  to  exist ; it  was  suggested  that  normal 
acceleration  might  also  serve  as  a piloting  cue,  although  it  would  not  likely 
be  so  important  as  pitch  attitude.  A final  assumption  was  that  the  pilot's 
primary  response  is  stick  force  F . The  resulting  pilot-vehicle  system  block 
diagram  is  shown  in  Figure  1;  the  locus  of  closed  loop  roots  of  the  system  is 
indicated  in  the  generic  sketch  for  the  case  of  no  significant  nonlinearities 
or  control  dynamics,  and  with  the  short-period  approximation  employed  for 
airframe  dynamics. 


Pilot 

Airframe 

■ • 0 f 

x V 1 „ , 

F„ 

* 

y — Vp(i-) 

Yp(l»>  * KP 

a a 

(The  "Synchronous' 

' Pilot) 

S_e\ 

H • • | 

K <»->  - < 

<* 

<j->  ii“’  ■ : > 

' s 'o  Jw(  ( j- 

‘ ) ' * : sp-sp 

(jw)  - control  'feel  dynamic 


(..**)  - static  gain  of  * it-) 

r*  n s 


^9P,W*p) 


Figure  i. 


! \ From 

High  Frequency  Asymptote:  oA  - ^ J * 2 

Simplified  Pilot-Vehicle  System  Model 
for  Fully-Developed  PIO 


From  Reference  4 and  5 
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The  non- equalized,  no  time  delay  model  for  pilot  dynamics — the  so-called 
"synchronous  pilot" — has  since  been  criticized  by  other  researchers  as  being 
too  unsophisticated.  One  must  be  careful  not  to  confuse  "unsophisticated" 
with  "elegant."  The  synchronous  pilot  model  is,  in  this  author's  opinion, 
the  elegant  product  of  very  sophisticated  reasoning.  It  was  developed  by  the 
same  people  (and  organization)  responsible  for  the  development  of  the  engi- 
neering theory  for  pilot  dynamics  (Ref.  b) ; the  pure  gain  form  of  the  model 
was  believed  to  be  in  reasonable  agreement  witli  the  very  limited  time  response 
data  available  for  actual  P10  experiences.  The  model  contains  two  features 
that  are  probably  essential  to  the  understanding  of  PIO: 

1.  The  overwhelming  importance  to  the  development  of  PIO  of  pitch 
attitude  loop  dynamics  is  explicit  in  the  model. 

2.  The  frequency  at  which  PIO  is  likely  to  occur  is  largelv  determined 
by  the  pitch  attitude  loop  closure. 


Each  of  these  points  will  be  discussed  later  in  this  report. 


Ashkenas  (Ref.  5)  used  the  linear,  pilot-vehicle  representation  of 
Figure  1 to  make  his  point  that  PIO  cannot  occur  in  the  absence  of  control 
system  dynamics--l i near  or  nonlinear;  that  is,  the  system  of  Figure  1 can 
never  become  unstable  for  real  aircraft  (witli  the  possible  exception  of 
closely-coupled  canards).  However,  his  basic  assumption  was  that  PIO  results 
solely  from  pitch  attitude  tracking;  lie  further  implied  that  for  PIO  to  occur 
it  is  necessary  that  either  a locus  crossing  of  the  imaginary  axis  must  occur 
in  Figure  1 as  a result  of  higher  order  system  dynamics  or  a limit  cycle  in 
the  control  of  pitch  attitude  must  occur  as  a result  of  control  system  non- 
linearities.  In  either  case,  Ashkenas'  theory  underscored  the  importance  to 

PIO  mechanics  of  1/T-  — a parameter  usually  ignored  in  experimental  work  of 

02 

that  period  and  known  to  be  of  major  significance  to  the  pitch  tracking  task 
handling  qualities.  This  viewpoint  was  the  genesis  of  the  PIO  criterion, 
based  entirely  on  pitch  control,  proposed  in  B22  as  a specification  for 
aircraft  design. 
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A'Harrah  (Ref.  7)  disagreed  sharply  witli  Ashkenas  on  the  significance  to 

I'll)  of  l/'l'o  , the  pitch  attitude  cue,  and  control  system  dynamics  or  non- 

linearities.  lie  contended,  based  on  experimental  results  obtained  using  a 

"g-seat"  (visual  plus  normal  acceleration  cues)  in  which  u>  , 0 and  F /g 

sp  sp  s 

were  varied  at  constant  l/Tp  , that  the  dominant  cue  is  normal  acceleration — 
not  pitch  attitude.  The  ensuing  argument  was  never  entirely  resolved. 
However,  Cornell  Aeronautical  Laboratory  flight  tests  (Ref.  17)  have  since 
established  that  P10  does  not  necessarily  result  from  short-period  dynamics 
which  force  to  lie  in  the  unstable-half  plane  (see  Figure  1).  In  addi- 
tion, as  Ashkenas  noted,  the  violation  of  A'llarrah's  PIO  boundaries  is  not 
suit  icient  for  1M0  to  he  created.  We  might  suspect,  therefore,  that  neither 
viewpoint  is  entirely  correct. 

Recent  flight  tests  (Ret.  8)  encountered  a serious  PIO  in  landing  flare 
with  an  all-linear  simulation  of  the  YF-17.  This  is  hard  evidence  that 
control  system  nonl inear i t ies  may  not  be  necessary  to  t lie  development  of  PIO. 
Control  system  phase  lags  were,  however,  identified  in  Reference  8 as  the  PIO 
cat  a 1 yst . 

There  are  two  obvious  difficulties  with  the  synchronous  pilot  model: 

1.  It  provides  no  insight  into  the  mechanisms  by  which  a fully 
developed  PIO  can  be  initiated. 

2.  The  assumption  that  the  pilot's  time  delay  in  single  loop  tracking 
can  approach  zero  may  be  indefensible  on  philosophical  grounds. 

The  first  deficiency  limits  our  ability  to  predict  PIO,  to  understand  it,  or 
to  develop  design  rules  for  avoiding  it.  The  second  bears  comment.  In 
tracking  of  periodic  Inputs,  it  has  been  observed  that  the  pilot's  time  delay 
appears  to  approach  zero  with  practice  (Ref.  9);  this  has  been  attributed  to 
the  pilot's  adaptation  of  a more  sophisticated  form  of  control  called  pre- 
cognitlve  (knowledge  of  input)  tracking.  This  conclusion,  however,  implies 
that  the  pilot  knows  what  his  time  delay  is;  otherwise,  he  has  no  wav  of 
eliminating  it.  This  is  assumed  here  to  be  a philosophical  impossibility. 

Are  you,  for  example,  reading  this  passage  now  or  did  you  read  it  yesterday 
and  only  believe  you  are  reading  it  now? 
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The  assumption  of  synchronous  control  behavior  is  based  on  the  assumption 
that  the  human  pilot,  in  single  loop  tracking,  may  be  represented  as  a single 
input/single  output  device  as  shown  in  Figure  1.  He  can,  of  course  (Ref.  6). 
Problems  arise,  however,  when  we  attempt  to  interpret  physical  data  against 
the  theoretical  framework  provided  by  such  a model.  Reference  10  challenges 
the  validity  of  the  servo  model  for  pilot  dynamics  (Ref.  6)  as  a description 
for  the  physiological  processes  involved  in  human  pilot  dynamics.  It  is 
suggested  in  Reference  10  that  when  the  pilot  is  modeled  as  a multi-sensor 
device  a unity  of  understanding  will  emerge  to  link  physical  response  data 
with  a rational,  theoretical  description  for  pilot  dynamics.  It  is  suggested 
that  the  resulting  model  will  be  completely  consistent  with  the  servo  model 
of  Reference  6 in  all  phenomenological  respects. 

A consequence  of  the  theory  of  Reference  10  is  that  the  input  to  the 

human  pilot  is  most  appropriately  represented  as  an  information  vector  (in  a 

manner  completely  consistent  with  modern  state  space  analyses  of  system 

dynamics).  Thus,  when  pitch  attitude  is  displayed  to  the  pilot  (in  the 

absence  of  motion  cues),  pitch  attitude  and  pitch  rate  should  be  used  as 

inputs  to  the  pilot  model  if  a mathematical  simulation  of  attitude  control  is 

to  be  performed  (Figure  2).  The  transfer  devices  labeled  G (ju>)  and  G (jw) 

q p 

are  intended  to  represent  the  dynamics  involved  in  the  visual  detection  and 
control  of  q(t)  and  0(t),  respectively.  These  elements,  especially  G^(jw), 
may  be  nonlinear;  however,  as  a first  approximation  they  may  both  be  replaced 
with  constant  gains  (Ref.  10).  When  the  q-loop  is  closed  (analytically),  the 
resulting  model  for  F^/0^  is  equal  to  Y^(jw)  as  given  by  the  servo  theory 
(Ref.  6). 


Pilot  Dynamic  F laments! 


•Ge(Ju»)  - ? 

(q^(  )u>)  - neuromuscular  system  dynamics 

Figure  2.  Multiple  Loop  Pilot-Vehicle  System  Model  for 
Pitcli  Attitude  Control  (no  Motion  Cues) 


Tin*  model  ol  figure  2,  when  compared  with  the  conventional  servo  model 
ot  Reference  h,  may  seem  needlessly  complex.  It  Is  not.  In  fact,  for  reasons 
discussed  In  Reference  10.  The  purpose  of  the  model  was  not  to  contribute  to 
the  advancement  ot  the  computational  aspects  ot  man/machlnc  systems  analyses, 
hut  to  enable  the  forecasting  or  explanation  ol  the  qualitative  relationships 
found  in  tasks  Important  to  vehicle  handling  qualities.  As  an  example, 
consider  how  the  model  ol  figure  2 might  he  used  to  devise  an  explanation  for 
the  physical  observation  that  the  time  delay  between  1'  (t)  and  0 (t)  approaches 

NO 

zero  in  sine  wave  tracking  experiments.  Consider  the  following  tails: 

1.  With  a periodic  input  to  a stable  servo  system,  the  average  system 
error  will  be  zero  when  the  controller  is  a pure  gain. 

2.  q (,  [ ) leads  0 til  bv  '*0  degrees;  the  correspond  1 ng  t line  advance  is 

e e 

n/2,0  seconds  (.ntg  frequency  ol  the  periodic  input), 
k 

1.  Control  ol  only  q(t)  ensures  control  ol  t'(l)  except  lor  possible 
low  I requeue v drift  In  error. 


We  might  speculate,  therefore,  that  for  reasons  that  are  as  yet  unclear,  t lie 
human  pilot  really  only  controls  error  rate  fin  a sine  wave  tracking  experi- 
ment) in  the  cent  inuous,  closed  loop  sense.  He  can  null  drift  in  orrot  with 
control  amp  1 1 1 tide  modulation.  further,  it  his  time  delay  in  rale  tracking  is 

i 

’’  / -’o  k seconds  then  it  will  appear  that  he  controls  system  error  with  zero 
l I me  delay  - the  synchronous  pilot  --  II  one  insists  that  the  system  be 
structured  as  shown  in  figure  1.  This  is  an  out  irelv  plausible,  physical  and 
(heoiettr.il  explanation  I oi  synchronous,  precoguitive  tracking  with  visual 
inputs  t o the  human  controller. 

I'he  loregoing  remarks  arc  presented  I o suggest  that  there  is  an  allot 
native  approach  to  visualization  and  modeling  ol  human  pile!  dynamics,  and 
that  this  mav  he  usetul  lot  decoding  I’lO.  The  model  ot  figure  is,  by 
dot  tuition,  quantitatively  identical  I o the  servo  model  ol  Reference  h lot 
visual  tracking,  ol  pitch  attitude.  A in.  i |or  advantage  ol  the  model  I ot  the 
study  ol  Pl<>  is  the  maimer  in  which  il  mav  be  extended  t o account  lot  tin' 
cl  lists  ol  tin  > t ion  cues  (lineal  and  angulat  accelerations).  The  rules  lot 
this  are : 


w 


Kach  additional  cuo  forms  a path  parallel  with  those  lor  q and  0 
(Figure  2).  The  associated  dynamics  are  those  of  the  motion  sensor 
mechanism  with  an  input-adaptive  gain. 

The  response  for  each  motion  cue  path  is  connected  witli  the  signal 
pathway  to  the  neuromuscular  system  through  a "switch."  This 
switch  is  postulated  to  exist  within  a "central  processor."  Its 
nature  is  such  that  it  may  appear  to  permit  coordinated  control  of 
all  feedback  cues  or,  under  stressful  conditions,  permit  only 
mot  ion  cue  control. 


This  situation  Is  depicted  in  Figure  3 for  the  case  where  the  motion  cue  is 

assumed  to  he  normal  acceleration  ;\y  at  the  pilot's  location  (observe  the 

P 

sign  convention).  The  task  is  assumed  to  he  regulation  of  pitch  attitude  in 
the  presence  of  atmospheric  turbulence;  therefore,  no  command  input  need  he 
shown  on  the  figure.  Figure  3 suggests  that  if  tht  "switch"  were  set  such 
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Figure  ).  The  Multiple  Loop  Pilot  Model  for  Pitch  Attitude 
Control  with  Normal  Acceleration  Cue 


that  only  acceleration  is  controlled  in  a fully-developed  P10,  then  F would 

s 

appear  to  be  synchronous  with  0 when  a phase  lag  exists  in  0 (jw)  of  magni- 
tude 


,0°  + * 


-til 


(jwpio) 


The  equivalent  time  delay  is 


r 


a 


(jwpI0) 


(57.3  uipjcj) 


A typical  value  for  this  delay  would  be  about  0.1  to  0.15  seconds,  based  on 
available  PH)  case  history  data.  This  point  is  noted  as  a matter  of  interest 
and  to  establish  a timewise  consistency  between  the  model  of  Figure  3 and 
that  of  Ashkenas  (Figure  1). 


it. 


A COMMENT  ON  STICK  PUMP  INC  AND  MOTION  CUES 


in  the  remainder  of  this  report  it  will  be  assumed  that  pilot-felt 

normal  acceleration  a,  is  the  only  motion  cue  that  must  be  considered  in  the 

zp 

analysis  of  PIO.  The  assumption  that  in  a fully-developed  1' ltf  the  pilot 
primarily  seeks  to  control  az . ^ appears  to  be  consistent  with  available 
evidence,  although  hard  data  are  scarce.  The  T-38A  PIO  time  history  con- 
tained in  Reference  4 and  repeated  here  in  Figure  11  is  the  best  such  example 
available.  Time  history  data  for  the  YF-12  PIO  are  given  in  Reference  II; 
the  presence  of  control  system  nonlinearities  makes  the  interpretation  of 
these  data  very  dtfitcult,  however. 


There  is,  to  this  author's  knowledge,  not  a single  example  in  the 
aircraft  handling  qualities  literature  to  support  the  inclusion  of  rotational 
acceleration  0 as  a primary  piloting  cue  in  a closed  loop  VFR  tracking  task. 

0 was  postulated  in  Reference  2 to  be  of  signal  importance  to  aircraft 
handling  qualities  and  as  the  cue  most  closely  linked  with  the  control  stick 
pumping  phenomenon.  However,  stick  pumping  was  exhaustively  analyzed  in 
Reference  3 (unfortunately  never  released  for  publication)  where  it  was 
determined  that  no  case  for  0 as  the  pumping  cue  could  be  made  that  was 
consistent  with  flight  test  records  and  pilot  comments  for  carrier  approach 
and  landing. 


However,  it  was  observed  that  to  a good  approximation  the  amplitudes  of 
az^(t)  oscillations  for  the  F-4B  in  actual  carrier  approach  were  nearly 
constant  and  equal  to  a value  only  slightly  greater  than  the  perception 
threshold  for  spine-wise  linear  acceleration.  The  corresponding  0 amplitude 
was  about  ten  times  the  known  perception  threshold. 


The  explanation  offered  (Ref.  3)  for  the  pumping  phenomenon  in  carrier 
approach  was  that  the  pilot's  creation  and  monitoring  of  a small,  oscillatory 
a provided  him  with  an  adaptive  cue  for  the  rapid  and  sensitive  detection 
of  aircraft  settle  due  to  downwash  aft  of  the  carrier  ramp.  Without  this 
cue,  the  delay  time  in  the  detection  of  monotonic  a_  due  to  downwash  would 
be  prohibitively  large.  The  detection  of  settle  with  this  pilot-created, 
adaptive  cue  appears  to  result  in  a throttle  feedforward  command  to  counter 
the  tendency  to  depart  from  the  desired  flight  path.  Subjective  support  for 
this  conclusion  is  offered  in  the  pilot  commentary  of  Reference  12. 


The  importance  of  stick  pumping  to  the  PIO  phenomenon  is  that  it  appears 
to  help  unify  our  understanding  of  how  skilled  pilots  may  selectively  choose 
those  aircraft  responses  that  are  most  directly  related  to  their  concept  of 
the  flight  control  task. 

I(  will  be  assumed  without  further  comment  that  the  pilot's  response  is 
stick  force  F . The  merits  of  the  assumptions  regarding  az  and  F^  may  lie 


judged  by  the  resulting  success  of  the  PIO  theory,  contained  in  the  next 
chapter,  as  a framework  for  understanding  the  PIO  case  histories  discussed  in 
Section  VI . 


C.  A PILOT  MODEL  FOR  PIO  ANALYSIS 


The  structural  form  for  the  pilot-aircraft  system  dynamics  is  assumed  to 
be  the  same  as  that  shown  in  Figure  3b.  In  the  remainder  of  this  report  the 
neuromuscular  dynamics  will  be  neglected;  i.e., 

Om(Jw)  = 1 


A further  assumption,  consistent  with  the  theory  of  Reference  10,  is  that  the 
pilot's  acceleration  channel  dynamics  are 


!a«“>  " Ka(%e 


- r. 


i.e.,  an  input-dependent  gain  with  a time  delay.  The  resulting  PIO  model  is 
shown  >n  Figure  4.  The  visual  tracking  portion  of  the  pilot's  dynamics 
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Figure  4.  A Pilot-Aircraft  Model  for  PIO  Analysis 


Y (ju>)  is  chosen  to  be  consistent  with  the  conventional  servo  model  repre- 
sentation; it  is  assumed  that  the  parameters  of  Y (1«>)  can  be  selected  to  be 

P 

consistent  with  the  "adjustment  rules"  contained  in  Reference  6. 


The  functional  dependency  between  K and  a.,  was  inspired  bv  a model 

a zpe 

for  visual  detection  of  signal  coherence  introduced  in  Reference  9 as  a model 
for  precognitive  tracking  of  a periodic  input.  It  arises  from  the  assumption 
that  the  process  by  which  a human  identifies  a stimulus  is  similar  to  the 
operation  of  a linear,  optimal  filter.  Based  on  optimal  filtering  theorv,  a 
model  for  the  acceleration  portion  of  the  pilot  model  is 

«a(J“>  = = Ka0  ( ~i/ir2)  (c°^R*aH"ra>w 

Pe 


where  and  A are  the  center  frequency  and  "width"  of  the  az^(t)  power 
spectrum  as  indicated  in  the  sketch: 


dfc) 


When  a2p(t)  has  a broadband  spectrum,  A is  large,  A is  small  and  the  corres- 
ponding pilot  gain  F„/a_  * 0.  When  a,  (t)  has  a narrowband  spectrum,  A 

4Pe  P 

is  small,  A is  large  and  Fs/az  is  maximum.  In  a time-wise  sense,  a narrow- 

Pe 

band  signal  is  one  that  is  nearly  periodic — much  like  a sinusoid  having  a 
randomly  varying  amplitude  and  phase.  A closed-loop,  pilot-aircraft  system 
that  is  "resonant,"  or  that  has  a lightly-damped,  dominant  mode,  might  have 
pitch  attitude  and  normal  acceleration  responses  with  narrowband  spectra. 
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The  time  delay  1 ^ is  unknown;  it  may  be  a function  of  the  center  fre- 
quency i ^ or  it  may  be  a constant.  As  an  interim  model,  it  is  suggested  that 
xa  = 0.25  seconds;  this  is  the  value  used  in  all  of  the  PIO  examples  in 
Sect  ion  V ! . 


MODEL  PROPERTIES 


It  is  postulated  that  in  pre-PIO  flight,  the  pilot-vehicle  system's  loop 

dynamics  are  dominated  by  the  pitch  attitude  loop  closure  (0,q  •+  F^);  that 

is,  we  postulate  that  prior  to  PlO-initiation  either  (1)  the  power  spectral 

density  of  a.,  (t)  is  broadband;  therefore,  K (a,  ) -*  0 and  there  is  no 
P a zp 

substantial  acceleration  tracking,  or  (2)  the  mode  ''switch"  is  usually  set  to 
activate  the  pitch  attitude  loop  to  the  near-exclusion  of  acceleration 
tracking. 

As  was  stated  above,  the  visual  portion  of  the  pilot  model,  Y^(jo>), 
is  assumed  to  be  completely  consistent  with  the  servo  model  for  pilot  dynamics 
as  documented  in  Reference  6;  this  portion  of  the  model  may  be  refined  for 
applications  at  the  user's  discretion  if  lie  wishes  to  account  for  the  effects 
of  motion  (az^  and  0)  on  the  model's  parameters.  Rules  for  doing  this  have 
not  been  thoroughly  researched;  data  are  sparse  and  somewhat  contradictory. 
Generally,  it  seems  to  be  true  that  the  pilot  gain  for  pitch  tracking  in 
flight  is  slightly  less  than  that  measured  in  fixed-base  simulation;  the 
svstem  crossover  frequency  is  correspondingly  reduced.  The  matter  is  not 
considered  to  be  of  particular  importance  to  present  purposes  and  will  not  he 
discussed  further. 

The  properties  of  the  mode  "switch"  (Figures  3 and  4) — if  it  exists — are 
unknown.  It  should  be  noted  that  this  switch — hypothesized  to  simulate 
higher  processes  within  the  central  nervous  system — may  be  redundant  if  the 
validity  of  the  input-adaptive  acceleration  gain  is  accepted.  It  will  be 
retained  in  this  report  because  both  possibilities  should  be  admitted  until 
better  information  becomes  available  and  because  it  makes  no  tangible  dif- 
ference to  our  final  result.  The  function  of  the  switch  is  to  provide  a 
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rational  connection  between  system  and  pilot  dynamics  resulting  mainly  from 
pitch  attitude  tracking  and  those  due  malnlv  to  normal  acceleration  tracking. 
It  is  assumed  that,  because  of  pilot  tension  (i.e.,  attentiveness,  anxietv, 
etc.),  and  flight  control  task  constraints,  the  pilot  switches  from  pitch  to 
acceleration  control  In  a logical  manner;  it  is  further  assumed  that  this 
switching  can  be  modeled  as  Indicated  in  Figure  4 with  an  appropriate 
switching  logic. 


An  interesting  and  plausible  model  for  mode  switching  is  suggested  hv 
Reference  9;  the  authors  suggest  that  an  input  stimulus  to  a pilot-vehicle 
system  may  be  considered  to  be  "subjectively  predictable"  when  v,  the  index 
of  subjective  predictability,  - 0.3. 


It  we  assume  that  the  coherence  properties  of  a stimulus  are  svnt  lies  i ,-ed 
within  the  higher  centers  (i.e.,  not  by  the  stimulus  sensor),  then  it  is 
reasonable  to  assume  that  the  same  model  t or  subjective  predictability  should 
applv  to  both  motion  and  visual  cues.  Speculating  further,  one  may  expect 
that  in  a state  of  incipient  I'lO  the  pilot  will  begin  to  emphasize  a., 
control  because  it  is  "predictable"  and  under  normal  circumstances  would 
therefore  be  controllable.  Thus,  v - 0.1  might  serve  as  the  required  switching 
logic.  Note,  however,  that  v is  determined  from  the  properties  of  tin'  normal 
accelerat ion  spectrum  for  pitch  attitude  tracking,  only.  It  follows  that 
pilot-aircraft  system  conf igurat ions  that  are  prone  to  closed  loop  "reso- 
nance" in  pitch  attitude  tracking  tasks  are  automat ica 1 Iv  i’10  candidates. 


A tinal  point  about  ti>e  switching  hypothesis  should  be  considered.  When 
an  actual  1*10  is  initiated  in  a real  airplane,  the  pilot  hasn’t  the  luxury  el 
coutemplat ing  system  responses,  assessing  whether  he  is  in  a I'lO  mode,  and 
selecting  an  appropriate  control  action.  Probably,  many  pilots  will  think 
there  is  something  the  matter  with  the  airplane  fe.g.,  SAS  failure)  and  will 
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continue  doing  whatever  it  was  that  catalyzed  the  oscillations — and  this  is 
entirely  reasonable.  Titus,  If  the  airplane  is  felt  to  sink  uncontrollably 
following  i’10  initiation  in  flare,  the  pilot  will  pull  back-stick  to  arrest 
the  normal  acceleration.  In  terms  of  the  switching  hypothesis,  this  suggests 
that  when  danger  due  to  flight  path  departure  is  lodged  to  be  immediate,  the 
attitude  control  mode  will  probably  become  secondary;  it  is  conceivable  that 
in  a fuilv  developed  PIO  such  as  experienced  on  several  occasions  with  t he  T- 
)8A  (Ref.  1),  the  pilot  "switches  into"  control  of  normal  acceleration  and 
cannot  easily  sat  1st v his  internal  logic  for  switching  back  to  pitch  attitude 
control.  This  limiting-case  description  for  fully  developed  PIO  may  not  be 
strictly  correct;  it  will,  however,  almost  certainly  provide  a system  of 
necessarv  conditions  which  must  be  satisfied  it  a fully-developed  PIO  is  to 
be  possible. 


JO 
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SECTION  IV 

\ no  theory 


A.  IMO  CATEGORIES 

It  tho  Pandora's  box  ot  possible  P!0  causes  and  ot toots  t omul  in  tin1 
litoiaturo  is  t o ho  svstomat  ioal  lv  addvossod,  cataloged,  and  understood,  it 
is  invi’ssatv  th.it  wo  consldet  tho  mechanisms  bv  which  I'lO  is  oroatod.  Poi 
this  purpose,  two  broad  olassos  ot  PIO  ait'  dotiued  as  follows: 

o I'vpo  1 — PIO  Indueod  bv  dynamic  response  ot  tho  olosed  loop  pilot 
control  ot  pitch  attitude. 

o I'vpo  l l — PIO  induced  bv  non-tracking  control  or  disturbance. 

It  was  implied  in  Section  III  that  the  model  tot  pilot  aircraft  dynamics 

shown  in  Figure  -*  was  intended  tv'  address  Tvpe  1 PIO.  It  appears  taint  is 

assumed!  that  the  mode  l ntav  also  bo  applied  tv'  the  analysis  ot  I'vpo  11  PIO 

provided  tliat  tin1  "switch"  is  permitted  t o initiate  a.  control  starting 
1 P 

from  a no- track  ins  condition. 

The  control  or  disturbance  inputs  required  tv'  initiate  a I'vpo  11  PIO 
would  be  dosoribable  as  "abrupt”  and  ot  nui.cn  i t tide  sufficient  tv'  excite  anv 
resi'iiant  stick-tree  dynamic  modes  v't  the  aircraft  system. 

11.  TYPE  l PIO 


Phase  t:  Pitch  Attitude  Svstom  Pvnamics 


It  is  postulated  that  Type  1 PIO  begins  with  highly  resonant  closevl 
1 v'v'p  vlvnamics  result  1 ng  t rom  pitch  attitude  tracking,  only.  It  t ho  pilot- 
vehicle  system  can  be  shown  to  be  nonresonant,  then  I'vpo  l PIO  is  unlikely. 
The  relevant  system  block  diagram  is  that  v't  Figure  •*  with  tho  switch  sot  tv' 
preclude  a.  feedback  control.  The  pilot  model  tor  ft  control 
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is  parameterized  according  to  the  "adjustment  rules"  of  Reference  6 (modi- 
fied, it  desired,  to  account  tor  motion  cue  effects  on  closed  loop  properties), 

l'he  describing  function  parameters  K , , T,  , and  T.  and  the  closed  loop 

pel  1 

parameters  > (crossover  trequencv),  phase  margin  and  gain  margin  are  all 

V.’ 

dependent  upon  the  aircraft  dynamics  ■ /F  fj.);  the  process  hv  which  these 
are  selected  is  iterative  and  the  selection  process  is  not  necessarily 
unique.  For  present  purposes  ot  developing  and  demonstrating  a theory  for 
l’ 10,  we  need  not  he  too  ftissv  about  this  process;  for  an  analysis  of  possible 
I'lO  problems  with  a new  aircraft  for  which  no  flight  tests  have  been  made, 
thi'  parameter  i/at  ion  of  Y ( ju>)  is  extremely  important  and  must  be  done  with 
great  care. 


Following  parameter  selection  for  Y (i.)  the  closed  loop  dynamics  must 
be  computed: 


t) 

Ypljui)  - (Ju>)  , (|u>) 

a * Fs  tSe 

■ > . U = 6 , 

1 + Yp(Ju))  ~ (Ju)  ^ (ju>) 


where  it  is  understood  that  t)  (t)  represents  an  equivalent  command  input  to 

o 

the  system  of  Figure  4 due  to  (vertical)  w-gusts: 
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Vhe  gust  response  transfer  function  is  (employing  the  short-period  approxi- 
mat  ions) 
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The  turbulence  power  spectral  density  4>w  w,(co)  can  be  chosen  as  the  Dryden 

p p 

or  Von  Karman  or  any  other  model  that  may  best  represent  the  specific  flight 
condition  to  which  this  theory  Is  to  be  applied.  The  corresponding  0 

c 

power  spectral  density  is 


(j«>  *WgWg(u,) 


Following  closure  of  the  pitch  loop,  the  (uncontrolled)  acceleration  response 
power  spectral  density  (due  only  to  closed  loop  control  of  pitch  attitude) 
mav  be  obtained  as  follows: 
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A typical  pitch  attitude  loop  closure  will  yield  [ 0 / 0 0 ( j io ) j ' l (0  db)  for 
*•  u>c  ( the  crossover  frequency);  there  may  be  a closed  loop  amplitude  rise 
at  associated  with  small  closed  loop  damping  ratio.  At  frequencies  greater 
than  ujc  the  closed  loop  response  will  be  heavily  attentuated  except  (possibly) 
at  frequencies  corresponding  to  lightly  damped  modes  due  to  the  feel  system. 


actuator  dynamics,  structural  dynamics,  etc.  A typical  case  is 
in  the  following  sketch.  (The  zero  db  lines  are  intended  to  he 
live  for  the  case  where  a?  /0C  is  in  units  of  g/degree.) 


ill ust  rat  ed 
represent a- 


The  progression  I rom  pitch  attitude  loop  dynamics  to  the  normal 
acceleration  power  spectral  density  is  illustrated. 

These  sketches  Illustrate  a number  ol  points  that  are  ol  value  to 
understanding  t hi'  connection  between  pitch  attitude  dynamics 
potent  (a  I 1 or  Tvpo  I 1M0: 


and  t he 
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Any  feature  of  pitch  attitude  loop  dynamics  which  promotes  loop 
resonance  (amplitude  ratio  peaks  on  the  closed  loop  Bode)  is  a 
potential  cause  of  PIO.  This  would  include  poor  pitch  dynamics 
and  excessive  pilot  gain,  for  example. 


Values  of  1/Tq  that  approach  will  promote  peaks  in  the  az^ 
power  spectrum.  Small  1/Tg  relative  to  10  has  a slight  attenuating 
effect. 


Small  C will  promote  a„  resonance — particularly  when  u = w . 
sp  zp  r 1 sp  c 

Higher  frequency  modes  due  to  feel  or  control  system  dynamics,  for 
example,  are  heavily  attenuated  by  the  turbulence  power  spectral 
roll-off;  i.e. , they  are  not  easily  excited  by  turbulence. 

The  turbulence  spectrum  break  frequency  is  not  too  important  to 

the  development  of  a„  resonance;  small  w decreases  the  amplitudes 

zp  b 

of  az  (t)  but  has  minor  effect  on  subjective  predictability  of  az^ 
according  to  the  proposed  index. 

As  a rule  of  thumb  it  appears  that  large  1/Tg  and  small  closed 

2 

loop  dominant  mode  damping  ratio  promote  the  development  of  az^(t) 
with  narrowband  signal  qualities  and  will  therefore  promote  Type  I 


The  center  frequency  amplitude,  A,  and  width  parameter,  A,  may  be 
estimated  from  the  power  spectrum  of  az  by  any  convenient  means.  It  is 
suggested  that  a2  be  estimated  by  direct  integration: 
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Then  A = /2A  nnd  v = A/io^.  |Note  that  the  turbulence  power  spectral 

density  is  nKt  important  to  the  estimate  of  v provided  that  the  "shape"  of 
l'a;,u;/  (uj)  without  regard  for  amplitude  near  is  not  very  turbulence 
sens  it ive.  | 

If  't'azaz(u>)  *'as  no  discernible  center  frequency  for  all  realistic 
paramet erizations  of  Yp(Ju)  then  Type  I t'10  is  unlikely.  If  the  az^  spectrum 
has  more  than  one  candidate  "center  frequency,"  and  if  two  or  more  of  these 
resonances  satisfy  the  predictability  criterion  v - ().'3,  then  all  the 
correspond ing  center  frequencies  must  be  used  to  evaluate  the  phase  criterion 
to  in-  discussed  later  in  this  section. 

As  a practical  matter  it  may  be  impossible  by  analysis  to  certify  that 
tit*.'  pitch  attitude  system  cannot  be  made  to  yield  a suitablv  narrowband  a. 

-p 

response.  In  those  cases,  what  is  required  is  that  bounds  on  be  estab- 
lished (especially  its  maximum  value). 

A shortcut  to  the  above  analysis  methodology  is  possible;  it  is 
recommended  for  those  flight  control  tasks  where  uncertainties  exist  about 
t lie  importance  of  turbulence  to  the  basic  piloting  problem.  For  example, 
in  flare  t lie  importance  of  good  pitch  attitude  loop  dynamics  is  almost  an 
inarguable  point.  Hut  flare  is  Initiated  by  control  feedforward  from  the 
pilot;  there  may  be  no  substantial  turbulence.  We  may  assume,  however,  that 
pitch  regulation  will  be  necessary  to  compensate  for  errors  generated  by 
control  resulting  from  pitch  command  errors.  The  corresponding  pitch  command 
error  spectrum  is  unknown;  a conservative  estimate  would  be  to  model  it  as 
broadband  noise,  set  $0  t)  (u>)  1 in  the  above  equation,  and  estimate  v as 

before.  The  resulting  'I'aza7  (ui)  for  the  problem  previously  sketched  would 
be  as  shown: 
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Tills  sketch  suggests  that  the  estimation  of  A,  w and  A may  not  he  unique. 

K 

An  even  simpler  criterion  for  subjective  predictability  could  he  based 
on  the  dominant  mode  damping  ratio  t for  closed  loop  pitch  attitude  con- 
trol. Until  better  data  can  be  developed  it  is  suggested  that  az  (t)  is 
subjectively  predictable  when 


The  corresponding  resonant  frequency  should  be  set  equal  to  the  dominant 

mode's  undamped  natural  frequency.  The  validity  of  this  criterion  should  be 

assessed  on  a case-by-case  basis.  Clearly,  there  may  be  combinations  of  ^ , 

1 /In  • c.  . and  u)  where  it  will  be  inadequate. 

*sp  sp 

If  control  or  feel  system  nonlinearities  exist  that  have  significant 
effect  on  closed  loop  pitch  dynamics,  then  the  estimation  procedure  for  v and 
kJ  must  be  modified.  No  attempt  will  be  made  here  to  explain  how  this  mav  be 
done.  Describing  function  techniques  (Reference  11)  may  be  employed  or  the 
pilot-vehicle  system  may  be  simulated  on  an  analog  or  digital  computer  and 
an  "experimental"  estimate  made.  If  nonlinearities  have  a significant 
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effect  on  pitch  tracking  dynamics  then  (almost  by  definition)  the  handling 
qualities  will  suffer;  this  situation  might  require  correction  as  a prelude 
to  future  exploration  of  P10  behavior. 


If  severe  nonlinearities  exist  in  the  pitch  attitude  control  loop  that 
cannot  be  ignored  in  estimation  of  system  dynamics,  then  it  is  conceivable 
that  pitch  attitude  limit  cycles  will  occur.  If  this  does,  indeed,  happen 
the  limit  cycle  frequency  should  be  set  equal  to  a7^(t)  should  be  assumed 
to  be  subjectively  predictable,  and  the  analysis  for  possible  P10  should 
proceed  to  phase  2 (a7^  "*  dynamics).  The  presence  of  a limit  cycle  in 
control  of  d does  not  imply  that  PIO  will  exist,  by  the  PIO  definition;  it 
would  only  mean  that  the  first  necessary  condition  for  PIO  had  been  satis- 
fied. The  distinction  is  important. 

By  now  it  is  probably  clear  that  the  analysis  of  pitch  dynamics  is 
somewhat  artistic.  A suitable  pilot  model  Y^(jw)  must  be  parameterized  to 
represent  both  the  normal  and  possible  bizarre  forms  of  pilot  dynamics  in 
the  pre-PIO  phase  of  attitude  control;  this  is  not  exactly  a problem  amen- 
able to  state-of-the-ar t solutions.  A turbulence  model  must  be  selected  to 
realistically  simulate  the  frequency  nature  of  expected  gust  disturbances; 
the  break  frequency  is  the  most  important  turbulence  parameter  for  this 
analysis.  Fortunately,  it  appears  that  so  long  as  i'w>w  (w)  is  "represen- 
tative"  that  it  Is  not  too  critical  to  an  evaluation  of  PIO  potential. 
Finally,  the  power  spectral  density  predicted  for  az^  must  be  examined  to 
diagnose  whether  it  is  suitably  narrowband — and  therefore  subjectively 
predictable;  the  rules  for  doing  this  are  primitive,  at  best.  Thus,  the 
problem  of  modeling  the  pitch  attitude  system  is  very  much  one  where  engi- 
neering judgment  Is  indispensable  to  success.  Fortunately  for  some  of  us, 
judgment  is  not  quite  so  vital  in  phase  2 of  Type  1 PIO  assessment. 


Phase  2:  Normal  Acceleration  System  Dynamics 

I l the  closed  loop  control  of  pitch  attitude  is  found  to  produce  a 
subjectively  predictable  normal  acceleration  response,  then  it  is  postulated 
that  the  pilot  will,  at  some  time,  attempt  to  track  az^.  This  process  has 
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been  modeled  (Figure  4)  with  a switching  function  such  that  either  0 or  az  ^ 
is  controlled,  but  not  both.  It  was  previously  mentioned  that  the  hypothesized 
a_  -dependent  pilot  gain  K might  be  used  as  a switch  to  "turn  on"  tin* 

l\ 

feedback  control  of  a,  . 


It  is  further  postulated  that  when  the  pilot  begins  to  track  az^  his 

"rule"  for  system  closure  is  that  the  crossover  frequency  of  the  _a_z  loop 

must  equal  — the  resonant  frequency  of  the  az  (t)  response  due  to  pitch 

attitude  control,  alone.  It  is  assumed  that  he  selects  K to  establish  this 

(1 

condition.  The  appropriate  pilot  model  for  az^  control  is  simply 


It  is  suggested  that,  until  better  information  becomes  available,  an  appro 
priate  delay  is 


— 

T =0.25  seconds 
a 

This  value  appears  to  be  consistent  with  the  T-38A  P10  tracking  data  (Rets.  I 
and  4);  more  important,  however,  is  the  observation  that  it  seems  to  correlate 
the  results  of  the  numerical  studies  (Section  VI)  with  actual  PIO  experiences 
It  is  conceivable  that  t is  functionally  dependent  upon  A possibility 

is  that 


n 

T 5*  - • 

a 

This  dependency  arises  from  consideration  of  the  pilot's  motion  sensor  as  an 
optimal,  linear  filter.  This  model  for  Tfl  will  not  be  used  further  in  this 
report . 

A necessary  condition  for  the  existence  of  Type  I P10  may  be  simply 
stated  as  follows: 


’ 
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where  ui^  is  the  minimum  frequency  for  which 


LUm)  = ' N 


It  is  understood  that  L ( 1 to ) represents  the  linear  portion  of  the  a.,  *-  F 

zp  s 

loop  and  that  N is  a sinusoidal  describing  function  representation  of  the 

nonlinear  part.  That  is,  the  a.,  >•  F system  is  assumed  to  have  been 

/P  s 

arranged  in  the  following  form: 


The  Equivalent  az^  * Fs  System 

[When  nonlinearities  are  imbedded  in  inner  loops  within  the  feel  or  auto- 
matic control  systems,  this  formulation  can  be  difficult  to  accomplish.]  I 
no  important  nonlinearities  exist,  then 


N 1 


L(m)  = -5—  (ju) 
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In  general,  the  describing  function  N will  be  a function  of  oscillation 
frequency  and  amplitude;  i.e.,  N = N(m,  a2  ).  In  the  examples  contained  in 
this  report,  however,  N will  be  amplitude  dependent  only;  i.e.,  N - N(a.,  1 . 
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Observe  that  the  criterion  5-  is  merely  necessary  for  P10  and  is 
not  a guarantee  that  PIO  will  occur.  Further  discussion  of  this  point  will 
be  deferred  until  Phase  3 of  the  Type  I PIO  assessment  is  introduced. 

With  a linear  system  the  criterion  is  equivalent  to  specifying 

that  the  phase  margin  must  be  zero  or  negative  if  Type  I PIO  is  to  occur. 


The  interaction  between  the  0 -*■  F and  a,  -*■  F loops  is  illustrated  in 

s *p  s 

Figure  5 for  a lightly  damped  linear  system  (representative  feel  system  and 

actuator  dynamics  are  shown;  a first  order  Pade  approximation  is  assumed  for 

the  time  delay).  Figure  5 is  intended  to  be  a generic  sketch,  only,  and  is 

not  drawn  to  scale.  The  figure  suggests  that  when  the  "switch"  of  pilot 

control  from  0 to  a,  occurs,  the  a_  •*  F system  will  be  stable  when  the 
zp  zp  s 

resonant  frequency  of  0 -*■  F^  is  less  than  the  frequency  for  which  the  az^  -*■  F^ 
locus  first  crosses  the  imaginary  axis.  When,  however,  the  resonant  fre- 
quency of  0 -*■  F^  is  greater  than  the  value  at  which  the  az^  ■+  F locus 
crossing  occurs  then  a^  ■+  Fr  will  be  unstable  if  the  pilot  switches  from 

control  of  0 to  control  of  ar  at  that  condition.  This  is  the  linearized 

P 

description  for  PIO  initiation. 


Solid  curves: 

locus  of  closed  loop 
roots  for  9 * FH 

Dashed  curves: 
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Figure  5.  Closed  Loop  Dynamics;  0 - Fs  and  az , > Fs 
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The  contributions  of  the  feel  system,  control  system,  pilot,  and  air- 
f r.nne  dynamics  to  P10,  for  a linear  system,  are  conveniently  displayed  on  a 
plot  of  total,  open  loop  system  phase  angle  t versus  frequency: 


$(.|u>)  " (pilot  phase  lag)  f (feel  system  phase)  + 

(control  system  phase)  + (airframe  phase) 

<S  _ az 

'(’(Ju))  “ -57.Jtau)  + i (Jio)  + ) (ju>)~  degrees 

A typical  plot  of  d*  ( J‘o ) might  appear  as  shown  below  in  figure  b. 
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flgure  b.  System  i'hase  Components  for  a-.,  ‘ fs 


It  the  aZj,  resonant  frequency  resulting  from  0 k Ks  closure  is  ^ , then 
the  a.-p  * Fs  system  phase  margin  dmi  - 1H0°  + d>  ( JtoR  ^ ) Is  positive,  the  a;.p 
system  is  stable,  and  no  f 10  results.  for  toR  - ,oR , , the  sketch  illustrates 
a negative  phase  margin  with  a resultant  closed  loop  instability;  this,  by  the 
fit'  definition,  is  a P10.  II  it  can  be  shown  that  u'g  > uq  is  impossible, 
then  Tvpe  1 1’  10  cannot  occur  for  the  associated  flight  condition.  This 
phase  margin  criterion  is  oversimplified  it  one  wishes  to  be  rigorous. 
Nvquist’s  criterion  can  be  checked  to  determine  stability  of  the  a _•  loop, 
given  the  crossover  t requencv  (and,  therefore,  the  loop  gain).  In  general, 
this  appears  to  be  an  unnecessary  comp  1 i cat  ton;  it  ,.'R  is  greater  than  the 
minimum  t requencv  tor  which  it  then  a*  k Vs  is  unstable  it  closure  is 

made  at  frequency  n'R. 


• Fg  dynamics,  however,  we  cannot  prescribe  a priori  what  the  character 


When  feel  and  control  system  nonlinearities  are  known  to  exist,  their 
effects  on  az^  ->  system  dynamics  must  be  considered.  Their  effect  is  to 
introduce  amplitude-dependent  and  frequency-dependent  gains  and  phases  into 
the  control  loops  for  both  0 and  aZp-  It  was  already  noted  that  the  effect 
of  nonlinearities  on  9 -*•  F^  dynamics  should  not  generally  be  too  significant 
for  small  motion  amplitude  unless  the  basic  handling  qualities  are  poor.  For 
a2 

of  the  system  dynamics  might  be;  allowance  must  be  made  for  the  possibility 
that  surface  or  rate  limits  can  be  encountered  due  to  large  amplitude  excur- 
sions, for  example. 

Generally,  reasonable  sinusoidal  describing  function  models  can  be 
obtained  for  each  important  nonlinearity  and  the  system  dynamics  can  then  be 
formulated  in  the  series  model  shown  previously,  combining  all  the  linear  ele- 
ments into  L(jui)  and  representing  the  sinusoidal  describing  functions  by  N. 

Then  limit  cycles  of  the  az  -»■  F loop  are  easily  found  from  plotting  ampli- 

P s 

tude  versus  phase  angle  for  both  L(jw)  and  the  negative  inverse  of  N on  the 
same  diagram.  Points  of  intersection  between  these  two  curves  correspond  to 
frequencies  and  amplitudes  of  motion  at  which  L(jw)  = -1/N.  The  smallest  of 
these  intersection  frequencies  is  defined  to  be  i^.  A typical  case  is 
illustrated  in  Figure  7 for  the  case  where  the  feel  system  contains  an 
amplitude-dependent  nonlinearity  and  no  others. 

In  Figure  7,  the  oscillation  frequency  U)  is  the  parameter  of  the  gain- 
phase  plot  of  L(jw);  oscillation  amplitude  az^  or  pilot  control  amplitude 
parametizes  the  -1/N  plot.  Three  points  of  intersection  between  these  two 
curves  are  indicated.  Points  1 and  3 are  stable  limit  cycles;  point  2 is  an 
unstable  limit  cycle. 

Limit  cycles  found  from  a describing  function  analysis  such  as  that 
indicated  above  do  not  necessarily  represent  PIO  conditions . For  a PIO  to 
actually  occur  it  is  necessary  (but  not  sufficient)  that  the  resonance  fre- 
quency jr  of  the  pitch  closed  loop  be  equal  to  or  greater  than  that  of  the 
lowest  frequency  limit  cycle  . Actually,  a better  prescription  for 
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Figure  7.  Determination  of  Limit  Cycles  in  aZp  Fs 

nonlinear  systems  might  be  «>  - • < , where  t is  a "small"  number:  that  is 

L K 

a limit  cycle  can  be  initiated  if  the  system  is  initially  oscillating  at  a 
frequency  "near"  the  limit  cycle  frequency.  It  may  not  be  possible  to 
specify  what  a suitable  "nearness"  criterion  would  be  and,  as  a practical 
issue,  it  is  not  important  here  (to  do  so  might  require  the  development  of  a 
criterion  based  on  Lyapunov's  second  method,  for  example).  This  nuance  will 
be  ignored. 


If  > u).  then  a PIO  mav  exist  with  initial  frequency  u> ,, . The  oscilla- 
t ions  may  orbit  into  a stable  limit  cycle  or  continue  to  increase  in  either 
amplitude  or  frequency  (or  both).  It  will  often  happen  that,  for  the  linear 
approximation  to  the  nonlinear  system  dynamics,  the  frequency  for  = -180” 
will  be  . In  those  cases,  instability  of  the  linear  system  will  imply 
that  Pit)  exists  with  the  nonlinear  system;  the  use  of  the  describing  function 
analysis  is  redundant  in  such  cases  for  the  determination  of  PIO  unless  one 
wishes  to  estimate  the  PIO  stability,  frequency,  and  amplitude. 


Is 


Phase  3:  Acceleration  Amplitude 

The  final  necessary  condition  for  Type  I P10  is  postulated  to  be  that 
the  amplitude  of  az  oscillations  due  to  pitch  attitude  control  must  be 


greater  than  some  "critical"  value  (az  ) 


The  "critical"  value  (az  ) 


° zp  CR  “ — * — zp  CR 

is  assumed  to  correspond  to  a pre-set  pilot  tolerance  level  such  that  when 

the  amplitude  of  az  (t)  oscillations  is  less  than  (az  ) the  pilot  is 

p p L-K 

insensitive  or  indifferent  to  it. 


The  conventional  viewpoint  of  the  human  pilot  is  that  he  is  highly 

adaptive.  It  follows  that  (az  )_  may  be  task-dependent.  For  example, 

p LK 

oscillations  of  a,  that  are  tolerable  without  pilot  control  in  air-to-air 
ZP 

combat  may  be  totally  unacceptable  in  aerial  refueling  or  in  flare.  A 
logical  corollary  to  this  reasoning  is  that  (az  ) will  vary  among  pilots 
and  among  pilot  communities;  a student  pilot,  for  example,  might  be  hypo- 
thesized to  emphasize  az^  control  more  than  might  a service  or  test  pilot. 
Thus,  it  might  be  concluded  on  this  basis  that  the  student  is  more  suscep- 
tible to  PIO  than  an  experienced  pilot.  This  may  be  so,  of  course;  but  it 
need  not  be  so  for  the  stated  reason.  If  this  viewpoint  of  the  adaptive 
pilot  is  accepted,  then  a contradiction  must  be  admitted;  viz.,  an  explana- 
tion must  be  offered  for  why  less  experienced  pilots  do  not  catalyze  more 

PIO's  as  a result  of  a,  -*  F tracking  in  view  of  the  fact  that  the  a7  F 

zp  s zp  s 

loop  dynamics  are  almost  universally  poor  and  prone  to  instability.  Rather 
than  compound  this  potential  philosophical  felony,  the  writer  suggests  that 
perhaps  we  should  reconsider  the  adaptive  capabilities  of  the  pilot. 


It  is  proposed  that,  if  a_  (t)  is  of  magnitude  sufficient  to  be  con- 

P 

sciously  felt  by  the  pilot,  and  if  it  is  subjectively  predictable  by  him, 

then  he  may  attempt  to  control  it.  If  a realistic  in-flight  threshold  of 

a7  sensation  is  assumed  to  be  0.01  g (Ref.  14),  then  it  follows  that  (a? 

^-p  ^p  CR 

= 0.01  g for  any  flight  control  task. 
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Using  the  data  summarized  in  Reference  10,  it  is  possible  to  conclude 
that  as  a first  approximation  the  pilot  opinion  rating  of  a pitch  attitude 
tracking  task  is  given  by  the  formula 


3.  y 


~ Cooper-Harper  scale 


The  overall  validity  of  this  formula  and  the  theory  on  which  it  is  based  need 

not  concern  us  here.  If  the  formula  is  accepted,  then  it  follows  that 

aq  = 0.833  degree/second  is  the  maximum  acceptable  value  for  rms  pitch  rate 

(i.e.,  for  POR  = 3.5).  it  follows  that  will  be  both  detectable  bv  the 

Zp 

pilot  and  of  concern  to  him  when,  to  a first  approximation. 


(j^R)  o > 0.01 

I ri  H 


0.012 


whert  ]aZp/o  is  in  units  of  g/degree/seeond . This  criterion  can  only  be 
verified  with  flight  test  experience.  It  will  be  shown  in  Section  \'l  that 
the  validity  of  the  criterion  is  confirmed  by  the  numerical  examples  (which 
represent  the  major  existing  Pit)  data  base!). 


The  most  important  single  parameter  to  the  response  ratio  aZp/0  is 
shown  in  the  Appendix  to  be  1 / T . To  a good  approximation 


7—  (s)  » 


: . 1 . 2U0  i 

-tx  8 + ^ 8 + tv 


Thus,  the  present  PIO  theory  has  implicated  1/T.,  as  a central  parameter  in 
two  of  the  three  criteria  necessary  for  obtaining  a Type  I PIO  (the  first  and 
the  third).  (See  Appendix  for  information  concerning  the  approximations  above.) 


c.  type  II  PIO 


I 


’ 

The  tasks  addressed  here  are  those  in  which  closed  loop  control  of  pitch 
is  not  an  a priori  requirement  for  PIO  onset.  This  might  include  high-g 
maneuvers  of  an  open  loop  control  sort,  trim  malfunction,  system  transients 
resulting  from  SAS/CAS  start-up  or  shutdown,  etc.  The  common  thread,  how- 
ever, is  that  the  control  (or  disturbance)  must  be  "suddenly  applied"  and  of 
amplitude  sufficient  to  excite  the  stick-free  dynamic  modes  of  the  aircraft. 

The  analysis  requirements  for  the  investigation  of  Type  II  PIO  are  much 
simpler  than  those  of  the  last  section.  Basically  the  procedure  is  the  same  , 

except  that  there  is  no  messy  attitude  loop  closure  to  be  performed. 

it  is  assumed  that  the  potential  for  Type  II  PIO  may  be  determined  as 
f allows : 


1.  Compute  the  power  spectral  density  of  aZp  using  a normalized, 

broadband  noise  representation  for  Fs(jw) — to  simulate  the  required 
"abrupt"  character;  that  is 


<t’azaz(w)  = 


! a v 

s _£ 


t 


2.  Determine  the  resonance  frequency  and  the  subjective  predictability 

index  v from  a (w)  exactly  as  was  done  for  Type  I PIO. 
zz 


3.  If  v £ 0.3  then  Type  II  PIO  cannot  be  ruled  out  as  a possibility; 
then  it  is  necessary  to  determine  uil  for  L(ju>)  = -1/N(u>)  as  was 
done  for  the  acceleration  control  loop  for  Type  I PIO. 

4.  If  vj  0.3,  oiR  > wL  and  | aZp/P ( jwR) | >0.012  g/degree/second, 

then  Type.  II  PIO  should  be  considered  probable  for  the  given  flight 
condition  at  some  time  during  the  aircraft's  lifetime. 
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An  alternative  procedure  which  is  simpler  and  which  may  he  just  as 
satisfactory  is  to  determine  whether  any  stick-free  dynamic  mode  exists  which 
significantly  contributes  to  5a  az(“’)  ant*  has  a damping  ratio  - 0.2. 

As  a rule  of  thumb,  the  modal  frequency  must  be  less  than  about  10  rad/sec 
for  sucli  a mode  to  produce  significant  effect  on  $aza  (ui)  due  to  high  fre- 
quency attenuation  of  az.  If  such  a mode  exists,  then  the  az^  response 
should  be  considered  to  be  subjectively  predictable,  u>  should  be  set  equal 
to  the  modal  frequency,  and  the  additional  criteria  for  Type  11  P10  examined; 

viz.,  is  ■'  r-  w and  does  |az  / 0 ( j to  ) | exceed  0.012  g/degree/second?  If  so, 

K L p ' k 

then  Type  II  P 10  is  probable. 


The  simplified  predictability  criterion  of  C.  ; 0.2  was  selected  based 
on  the  T-18A  and  AAD-2  PIO  experiences.  Table  2 contains  data  extracted  from 
page  54  of  Reference  15.  The  airframe  and  feel  system  dynamics  are  shown  for 
the  flight  conditions  where  each  airplane  experienced  serious  PT0;  these  data 
are  given  for  the  original  and  the  modified  control  systems.  It  is  clear 


TABLE  2.  SYSTEM  DYNAMIC  DATA;  T-38A  AND  A4D-2 


Airp lane 

STI  TM  239-3 

T-38A,  aft  c.g. 

Douglas  LB-25452 
A4D-2,  c.g.  slight  aft 

Control  System 

Orig 

Mod 

Orig 

Mod 

Airframe  Dynamics: 

| wsp 

Stick  Fixed  < 

f ’ sp 

\ IJsp 

Stick  Free  t t 

( rsp 

7.0 

0.40 

9.8 

0.10 

7.0 

0.40 

7.5 

0.28 

7.94 

0.40 

8.8 

0.16 

7.94 

0.40 

4.8 

0.34 

Feel  System  Dynamics: 

(in  flight) 

(w  FS 

17.7 

24 

18.4 

19.2 

U FS 

0.23 

0.17 

0.22 

0.14 
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from  these  data  that  the  stick-free,  short-period  mode  probably  dominates  the 
character  of  az^(t)  responses  to  abrupt  control  or  turbulence  inputs.  The 
further  assumptions  are  made  that  both  the  T-38A  and  the  A4D-2  were  suscep- 
tible to  Type  II  PIO  (and  this  appears  to  be  a self-consistent  assumption) 
and  that  the  problems  were  cured  (or  at  least  alleviated)  by  the  control 
system  modifications  performed.  An  inspection  of  these  data  suggests  that 

o £,r  = 0.16  will  permit  Type  II  PIO  (A4D-2 , original  sitck-free 
damping  ratio). 

o = 0.28  is  sufficient  to  eliminate  (or  reduce)  Type  II  PIO 

tendencies  (T-38A,  modified  stick-free  damping  ratio). 

It  follows  that  between  0.16  and  0.28  will  represent  a susceptibility 
boundary. 

Substantiation  for  this  criterion  is  available  from  Reference  16. 

Bobweights  were  identified  as  a possible  contributing  cause  for  PIO  for  those 
cases  where  they  degrade  the  stick-free  short-period  dynamics.  The  overall 
importance  of  stick  force  per  steady-state  g,  minimum  dynamic  stick  force  per 
g,  and  short-period  damping  ratio  was  cited.  The  suggested  correlation  is 
shown  in  Figure  2 of  Reference  16,  p.  161  and  in  Figure  8 below;  the  PIO 
rating  scale  used  is  that  given  in  Figure  36.  These  data  were  obtained  from 
a flight  test  program  documented  in  Reference  17.  It  is,  to  this  writer, 
clear  from  this  data  plot  that  stick  force  per  g is  a tenuous  basis  for  a 
PIO  criterion.  The  data  show  only  three  PIO  or  near-PIO  cases;  one  of  these 
lies  in  the  "safe"  region  of  Fs/g.  One  no-PIO  case  lies  on  the  1.4  lb/g 
boundary.  There  are  four  cases  with  PIOR  = 3-3.5  lying  in  the  region  of 
<V8>min  > 3'°  lb/g.  If  one  assumes  that  a correlation  exists  between  PIOR 
and  damping  ratio  regardless  of  Fg/g  then  it  is  clear  from  these  data  that  a 
damping  ratio  of  about  0.2  should  represent  a PIO  boundary.  This  is  indicated 
in  the  figure.  This  damping  ratio  should  be  interpreted,  in  general,  as  that 
of  the  resonant  mode;  i.e.,  as  ?R.  When  the  variations  of  PIOR  with  Cr  shown 
above  are  interpreted  against  present  theory,  the  data  are  completely  consistent. 
These  data  potentially  represent  both  Type  I and  Type  II  PIO  conditions  due  to 
the  nature  of  the  flight  test  procedures  used;  the  flight  tests  did  appear 
to  emphasize  open-loop  maneuvering  to  a considerable  extent,  however.  Since 
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Short-Period  Damping  Ratio,  r.Sp 


Figure  8.  Effect  of  s ,,  and  Stick  Force  Level  on  PIO 
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has  been  postulated  to  he  a necessary,  hut  not  sufficient,  condition 
■ II  PIO,  (t  Is  reasonable  to  expect  that  the  three  non- PIO  cases  tested 
^ O.,!  did  not  satisfy  the  additional  Type 
| az  / 0(  | ui^)  | > 0.012).  Similarly,  the  one  data  point  with 
2 and  PI  OR  > 2.5  may  possibly  result  front  tendencies  toward  develop- 
a Type  I PIO. 


PIO  criteria  (a..  *■  F 

P 


It  should  be  noted  that  the  YF-12  and  the  YF-17  fas  simulated  bv 
OAl.SPAN)  both  experienced  serious  PIO.  The  short-period  damping  ratio  for 
the  YF-17  is  estimated  to  he  almost  0.1>,  while  that  ot  the  YF-12  is  about 
0.4  — both  in  the  PIO  coni'  igurat  ions . The  interpretation  of  these  data 
against  criteria  based  onlv  on  F/g  Is  impossible.  Both  cases,  however,  are 
readllv  understood  when  compared  against  the  PIO  criteria  suggested  a;;  a 
result  ot  the  I'vpe  I and  l'vpe  It  distinctions  lor  Pitt.  These  case  histories 
will  be  discussed  in  some  detail  in  Section  VI. 


SECTION  V 


SUMMARY  OK  RULES  KOK  Pit)  ASSESSMENT 


This  chapter  is  presented  as  .1  summary  and  compendium  of  rules  for  the 
assessment  ot  aircraft  and  flight  control  systems  with  respect  to  their 
potential  for  1'1()  encounter.  It  is  a condensation  of  the  theory  of  the  last 
chapter.  It  is  assumed  that  the  flight  conditions  and  aircraft-control 
system  configurations  to  be  investigated  are  known  and  that  all  the  dynamic 
models  required  for  application  of  the  theory  have  been  measured  or  estimated. 


A.  TYPE  I PIO  (INITIATED  BY  PITCH  ATTITUDE  CONTROL) 

1.  Select  an  appropriate  model  for  piloted  control  of  pitch  attitude 
(e.g..  Reference  b) . 

2.  Close  the  pitch  attitude  loop. 

t.  Compute  the  power  spectral  density  Ia.  I"')  of  normal  acceleration 
at  the  pilot's  location  due  to  control  of  pitch  attitude  (l.e.,  with  no 
closed  loop  control  of  acceleration).  Use  a representative  model  for  vertical 
t urhu lence. 


4.  From  J.i  a l60)  estimate  w (the  resonant  frequency) — il  one  exists, 
z z K 

It  none  exists  for  any  realistic  choice  of  pilot  dynamics  (allowing  for  gain 

changes  or  drop-out  of  equalization,  for  example,  due  to  pilot  tension, 

etc.)  then  type  I PIO  is  unlikely.  If  w exists,  then  estimate  the  sub- 

K 

jectlve  predictability  index  v.  If  v > 0.3  go  to  step  4a;  otherwise  con- 
t i nuo  with  step  5. 


4a.  Estimate  the  resonant  mode  damping  ratio  4 . If  C ^ 0.2  then 

K l\ 

conclude  that  Type  I PIO  is  unlikely;  the  evaluation  can  be  continued  at 
step  3 if  conservatism  Is  required  or  if  one  has  little  confidence  in  the 
pilot  model  parameterization.  If  - 0.2  go  to  step  r). 


5. 


It  there  are  no  teel  or  control  system  nonlinearities  present  that 
may  have  significant  effect  on  pi  lot -vehicle  system  dynamics  in  high-g, 
oscillatory  states  of  motion,  continue  with  step  5a;  otherwise  go  to  step  5b. 

5a.  (Linear  system  dynamics)  1’lot  the  total  open  loop  system  phase 
angle  Bode  >Kjm)  for  the  a^  *■  F loop  dynamics.  1 aj ) will  he  the  sum  of 
phase  angles  due  to  the  pilot,  the  feel  system  dynamics,  the  control  system 
(or  CAS)  dynamics,  and  tiie  airframe  dynamics  (including  any  SAS  feedback 
loops).  The  pilot  phase  should  be  assumed  to  result  entirely  from  a 0.25 
second  delay  (i.e.,  -57.  3*0.  25  ).  If  t lie  phase  margin  180°  + 4*  ( 1 uv, ) > 0 
then  Type  I I’ll)  is  unlikely  and  tile  analysis  lor  Type  t i’10  is  complete. 

If  180°  + f>(.ju>K)  s ()  then  Type  I P10  is  possible  and  the  analysis  should 
proceed  with  step  6.  If  180°  + i}>(juy)  's  greater  than,  but  approximately 
equal  to  zero,  then  this  may  suggest  that  the  airplane  will  exhibit  PIO 
tendencies,  although  a fully-developed  PIO  may  be  unlikely.  In  that  case, 
it  becomes  a matter  of  judgment  as  to  what  course  to  follow  in  a design  or 
development  program. 


5b.  (Nonlinear  system  dynamics)  On  a gain-phase  plot  of  l.(ju>)  and 
-l/N--the  lineitr  and  negative  inverse  of  the  nonlinear  portions  of  the 

* F^  system,  respec t ively — determine  the  smallest  frequency  uVj  at  which 


these  two  curves  intersect.  if  oj 


R 


Wj  , a Type 
U WR  < “i 


P10  is  possible  and  the 
analysis  must  continue  with  step  b.  If  w < then  Tvpe  1 Pill  is  unlikely 
and  the  analysis  for  Type  I P10  is  complete  (the  cautionary  note  of  5a 
applies  to  tills  inequality  determination,  however).  The  pilot  model  for  the 
analysis  of  az^  ■+  dynamics  should  be 


Fg 

a 2 


Pe 


(jrn)  = Ka  e 'a'W 


it  is  suggested  that  x = 0.25  seconds  is  a satisfactory  choice. 

b.  If  |az  /O(jwg)^  0.012  g/degree/sec  then  conclude  that  Type  l PIO 
is  unlikely.  If  this  ratio  is  > 0.012,  then  conclude  that  Type  1 PIO  is  a 
very  strong  possibility. 


. 

1 


B.  TYPE  II  PIO  (INITIATED  BY  MANEUVERING  CONTROL  OR  TURBULENCE) 


1.  Compute  the  power  spectral  density  of  az^  for  the  stick-free 
airplane  dynamics;  assume  that  the  airplane  is  excited  by  a wideband  noise 
with  power  spectral  density  = 1 (to  simulate  "abrupt"  inputs).  That  is, 
assume 


cl  ^ 

dZ  2 


(u>)  = 


aZr 


(jw) 


2.  Continue  the  analysis  exactly  as  described  above  for  Type  I PIO, 
starting  at  step  4 and  replacing  "Type  I"  with  "Type  II." 

3.  Simplified  alternative  procedure:  if  the  damping  ratio  of  the 

dominant,  resonant  mode  of  az  /F  is  < 0.2,  then  Type  II  PIO  is  possible. 

P s 

Call  this  damping  ratio  t and  continue  the  analysis  as  described  above  for 

R 

Type  I PIO,  starting  with  step  5 and  replacing  "Type  I"  with  "Type  II."  If 

the  damping  ratio  is  >0.2,  then  conclude  that  Type  II  PIO  is  unlikely.  For 

conservatism  one  could  define  w„  as  the  dominant  mode's  damped  frequency  and 

R 

proceed  to  step  5 of  the  Type  I PIO  analysis  above. 


SECTION  VI 
NUMERICAL  EXAMPLES 

A.  THE  YF-17  (APPROACH  AND  LANDING:  CAS-ON) 

The  reader  should  note  that  the  YF-17  never  experienced  a PIO  problem  to 
this  writer's  knowledge.  The  NT-33A  variable  stability  airplane,  configured 
to  simulate  the  YF-17,  did  encounter  serious  PIO  difficulties  in  flare 
(Ret  . 8) . Calspan  diagnosed  the  problem  as  excessive  control  system  phase 
t!  lag  due  to  the  CAS  mechanization.  They  modified  the  (simulated)  control 

|[  system  dynamics  to  reduce  the  lag  contribution  to  longitudinal  dynamics  and 

found  that  this  eliminated  the  PIO.  This  occurred  prior  to  the  first  flight 
ot  the  prototype  airplane.  Based  on  the  variable  stability  airplane  simula- 
tions the  YF-17  flight  control  system  was  similarly  modified  prior  to  the 
> first  flight.  Throughout  the  rest  of  this  report,  this  PIO  case  history  will 

I be  referred  to  as  the  YF-17  PIO  without  further  qualification.  The  YF-17 

case  is  particularly  interesting  because  the  airplane-control  system  dynamics 
is  simulated  were  entirely  linear . 

System  Dynamics 

It  is  unfortunate  that  a quality  description  of  the  YF-17  airframe  and 
control  system  dynamics  was  not  available  for  analysis  and  presentation 
within  tliis  report.  The  only  data  available  to  this  writer  were  those  con- 
tained in  Reference  8;  this  included  only  the  pitch  attitude  transfer  func- 
tion with  no  stability  derivatives.  In  order  to  reconstruct  the  corres- 
ponding normal  acceleration  transfer  function  some  assumptions  are  required 
as  follows: 

I . From  Reference  8, 

ff  y M<s„  (s  + .84) 

(s)  = ~ ^ 

s-  + 2(  . 89)  ( 1 . 98)s  + 1.98- 

, 

f . . 
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Note  that  the  acceleration  dynamics  estimated  here  for  the  YF-17  are  not 
exactly  those  that  were  simulated  in  the  reference  8 flight  tests.  The 
simulated  acceleration  numerator  was  that  of  the  basic  NT-33A  — in  the 
absence  of  a direct  lift  control  capability  or  the  use  of  model-following 
techniques  (which  were  apparently  not  used).  The  present  analysis,  using 
estimated  YF-17  dynamics,  will  establish  that  according  to  present  theory  the 
YF-17  with  the  original  control  system  was  indeed  PlO-prone.  The  question  of 
in-flight  simulator  fidelity  will  be  further  examined  in  Section  VII  where  it 
will  be  concluded  that  the  Calspan  simulation  was  qualitatively  correct. 


6.  To  a first  approximation,  assume  the  feel  and  control  dynamics  to  be 
simplifications  of  those  given  in  Reference  8: 


^ (s)  = ^ (s) 


5 po 


<s) ' (Hrl  * t5-  <s> 


(s) 

Fs 


K*  [s~  + 2(.44)(ll)s  + 11*  ] (s  + 2) (s  + 2,3) 
[s2  + 2(.7)(4)s  + 42 ] <s  + -9)(s  + 5) 


~ original 


e , , K'  (s  + 18)  (s  + 2)  (s  + 2.3) 


(s) 


Fs  (s  + 10)  (s  + . 9)  (s  + 5) 


~ modified 


K.'  is  unknown  (and  not  required  in  this  analysis)  , 


7.  Summary: 


0_  , \ K(2) (2 . 3) [ . 44 , 11] 

Fs  ~ (0) [.89,1.98] (5) [ .7,4] 


~ original 


0 , , K(2) (2. 3) (18) 

Fs  (S;  ' (0) [.89,1.98] (5) (10) 


~ mod if ied 
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f!E  (s)  , K[.08>5.04](2)(2.3)[.44.in 

Fs  [.89,1.98]  ( . 9)  (5)  [.7,4] 


'•'original 


azP  , . K[.08,5.04](2)(2.3)(18) 

Fs  (S)  ~ [,89,1.98](.9)(5)(10) 


~ modified 


where  the  short  hand  notation  employed  above  is: 


(1/T)  means  (s  + 1/T) 


[C,u]  means  s2  + 2£ius  + u>2 


The  units  are 


radians,  ft /sec2  and  pounds. 
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Type  II  P10  Assessment 

The  simplified  criterion  of  the  last  section  (step  3)  allows  us  to 
immediately  conclude  that  Type  II  PIO  in  this  configuration  is  extremely 
unlikely.  The  only  second  order  mode  appearing  in  the  acceleration  transfer 


function  is  the  short-period  mode.  But  C = 

sp 

subjective  predictability  criterion  C ? 0.2. 

K 


0.89  is  far  in  excess  of  the 


Type  I PIO  Assessment 

Figure  9 is  a Bode  plot  of  the  airplane's  pitch  attitude  dynamics 
0/F  (ju).  If  we  assume  that  the  crossover  frequency  will  lie  between  2-4 
rad/sec,  then  it  is  clear  that  the  aircraft  dynamics  are  roughly  of  the  form 
K/s'  in  this  region.  As  a rule,  dynamics  of  these  sort  will  lead  to  lightly 
damped  closed  loop  oscillations  and  degraded  pilot  opinion  ratings.  An 
inspection  of  the  data  base  of  Reference  6 and  a modicum  of  iteration  sug- 
gests that  a reasonable  (not  bizarre!)  model  for  pilot  dynamics  in  pitch 
tracking  would  be 

Y (jw)  = K (2.5ju>+l)  e'’385jU) 

A Bode  plot  of  the  open  loop  system  dynamics  Y^fO/F^)  is  also  shown  in 
Figure  9.  Figure  9 indicates  that  the  absolute  maximum  crossover  frequency 
with  this  Y ^ ( j to)  is  3.3  rad/sec.  Accordingly,  = 2.9  was  selected  and  is 
assumed  to  be  consistent  with  what  would  be  measured  in  actual  flight;  this 
yields  a small  phase  margin  (about  16°).  Obviously,  even  small  increases  in 
pilot  gain  will  rapidly  degrade  system  stability.  This  result  appears  to  be 
consistent  with  the  evaluation  pilots'  comments  about  the  poor  pitch  handling 
qualities  of  this  configuration  in  flight  tests  (Reference  8). 


The  corresponding  closed  loop  dynamics  0/0  (Juj)  are  shown  in  Figure  l) 
for  ni  = 2.9.  Obviously,  the  closed  loop  system  is  extremely  resonant  at 
this  condition.  It  is  evident  by  inspection  that  the  resonant  peak  of  0/0 
will  dominate  the  az  power  spectrum.  The  corresponding  damping  ratio  for 
this  mode  is  approximately  0.03.  Thus,  by  t lie  simplified  criterion  tor 


Figure  9.  YF-17  Fitch  Attitude  Dynamics 

subjective  predictability  of  the  last  chapter,  it  must  be  concluded  that 

Type  I PIO  cannot  be  ruled  out  on  the  basis  of  pitch  control  handling 

qualities.  The  resonance  frequency  = 3.0  rad/sec  for  the  given  Y (ju>). 

R p 

More  pilot  lead  and  higher  gain  would  increase  ^ somewhat. 

K 


By  the  assessment  rules,  the  analysis  must  now  proceed  to  an  investi- 


gation of  stability  of  the  az  loop  when  the  pilot's  gain  is  adjusted  to 

make  The  total  az^  -*•  F system  phase  j'(juj)  versus  frequency  is 

plotted  in  Figure  10  in  accordance  with  the  rules  of  the  PIO  theory.  The 


pilot  time  delav  was  assumed  to  be  0.25  seconds.  At  u>  = 3.0  we  see  that 
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4’=  -205°,  180  + 4>  = -25°  (the  system  phase  margin)  and  we  see 
acceleration  closed  loop  is  unstable.  Thus,  longitudinal  I’ll' 
initiated  provided  that  the  pilot  attempts  to  control  a z ^ . To 
this  is  possible,  proceed  to  step  6 of  the  assessment  rules. 


that  the 
can  be 
see  whether 


The  ratio  [a  /©(3.0.j)|  = 0.031  g/degree/sec . Thus,  by  present  theory 
we  would  be  justified  in  concluding  that  Type  I PIO  would  be  likely  with  this 
airplane  and  control  system. 


The  actual  normal  acceleration  dynamics  simulated  with  the  NT-33A  yield 
az  /i  (3.0.0  = 0.0213  g/degree/sec  (Section  VII).  This  is  about  twice  the 
criterion  value  of  0.012;  on  that  basis  it  can  be  concluded  that  errors  in 
the  simulation  of  a.,  motion  amplitude  were  probably  of  no  consequence. 


The  PIO  frequency  and  amplitude  obtained  with  the  NT-33A  simulation  are 
unpublished.  It  is  known  from  informal  communication  between  the  writer  and 
t'al span  staff  members  that  the  PIO  frequency  occurred  at  approximately  1/2  cps. 
It  may  therefore  be  concluded  that  this  analysis  (and,  as  a consequence,  the 
present  theory)  is  supported  by  the  flight  test  results. 


But  what  of  the  YF-17  with  the  modified  control  system?  If  the  PIO 
theorv  is  to  be  acceptable,  it  must  show  that  PIO  is  no  longer  a probability 
in  this  flight  condition.  This  is,  indeed,  the  case. 


A Bode  phase  plot  is  shown  in  Figure  10  for  the  two  pre-filter  dynamics 
(orginal  and  modified)  flight  tested  by  Calspan.  It  appears  that  the  modi- 
f ied  control  system  pre-filter  decreases  total  system  phase  lag  in  the  region 
of  probable  pitch  system  crossover  by  45-70  degrees.  The  resulting  airplane 
dynamics  will  therefore  be  more  like  K/s  than  the  original  K/s'  form.  The 
handling  qualities  are  bound  to  be  improved  accordingly.  No  detailed  analysis 
of  pitch  tracking  is  necessary,  however,  beyond  the  comment  that,  because  of 
the  improved  control  system  dynamics,  the  chances  are  substantially  reduced 
for  obtaining  resonant  closed  loop  dynamics. 
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Figure  10.  YF-17  Acceleration  Control  System  Dynamics 


Consider  the  a„  ->  F_  phase  angle  criterion.  In  Figure  10,  4> ( i U3 ) is 


plotted  with  the  modified  control  system.  There  are  no  intersections  of  the 
-180°  phase  condition  for  u < 12.5  rad/sec  provided  ta  is  less  than  about 
0.3  seconds.  Type  I PIO  would  be  extremely  unlikely  with  the  modified  system. 


It  is  therefore  concluded  that  the  theory  of  this  report  is  substanti- 
ated by  the  YF-17  in-flight  simulation.  The  reasons  offered  in  Reference  8 
for  the  PIO  problem  are  confirmed,  also,  by  these  these  theoretical  results. 
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B. 


THE  T-38A  (M  * 0.  HI  AT  SEA  LEVEL) 


The  T10  problems  encountered  with  this  airplane  are  well-known  and  have 
been  given  considerable  engineering  attention  (e.g..  References  1,  4,  and 
19).  The  Pit)  experiences  with  this  airplane  were  we  1 1 -documented , consider- 
able analyses  were  conducted  (and  documented)  which  Isolated  the  device- 
centered  origin  of  the  flight  control  problem,  and  significant  control  system 
changes  were  made  to  the  airplane  which  generally  appear  to  have  solved  the 
problems;  thus,  we  have  a fairly  complete  before-and-af ter  hlstorv  ot  this 
airplane’s  Pit)  behavior. 

The  i’ll)  time  history,  cited  in  the  above  references,  is  reproduced  in 
figure  11  tor  the  reader's  convenience. 


Figure  11.  T-38A  1M0  Time  Hlstorv 


I 

System  Dynamics — Or  i g Lna  1 Control  System 

dynamics  of  the  airplane  plus  flight  control  system  are  complicated 

A normal  acceleration  bobweight 
Control  system  friction 

The  bobweight  had  a very  destabilizing  influence  on  the  short-period  dynamics; 
tills  can  be  seen  from  Table  2 (p.  18).  There  was  sufficient  control  system 

friction  (break-out  force  = 2 pounds)  to  ensure  that  the  bobweight  feedback 
was  almost  completely  masked  for  az  less  than  about  1 g (the  bobweight  gain 
was  2 lb/g).  The  effect  was  to  qualitatively  change  the  airplane  dynamics 
according  to  whether  az  was  small  or  large.  For  small  az  the  airplane's 
* dynamics  were  best  represented  by  the  stick-fixed  condition;  for  large  az  the 

stick-free  dynamics  were  appropriate.  The  change  in  dynamics  probably  occurred 
very  suddenly  with  changes  in  acceleration  amplitude. 

The  structure  of  the  airframe  plus  flight  control  system  dynamics  (SAS-off) 
is  shown  in  Figure  12.  The  bobweight  is  modeled  as  a feedback  of  az^  to 
stick  force.  The  pitch  attitude  dynamics  0/Fg(ju>)  are  shown  in  Figure  13  for 
the  two  cases  of  no  bobweight  and  full  bobweight. 

Type  I PIO  Assessment  — T-38A 

The  most  difficult  task  in  analyzing  the  Type  I PIO  behavior  of  the  T-38A 
is  to  select  suitable  models  for  pilot  dynamics  in  control  of  pitch  attitude. 
This  will  be  done  for  three  cases  which  sufficiently  cover  the  spectrum  of 
pilot-vehicle  system  dynamics: 

o Fully-adapted  pilot;  small  az  (Kg  = 0) 

o Non-equalized  pilot;  small  az  (Kg  = 0) 

o Fully-adapted  pilot;  large  az^  (Kg  = 2) 


The 


by 
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Figure  13.  T-38A  Pitch  Attitude  Dynamics;  Effect  of  Bobweight 

Reference  4 suggests  that  prior  to  PIO  initiation,  a reasonable  model  for  the 
fully-adapted  pilot  in  control  of  pitch  attitude  is 


Kp  e~-2s 

VS>  = s 


3.2 


+ 1 


The  lag  time  constant  was  selected  to  cancel  the  0/F  (s)  numerator  time 

s 

constant;  the  resulting  open-loop  dynamics  of  the  0 ->  F^  loop  are  therefore 
approximately  equal  in  form  to  K/s  in  the  region  of  expected  crossover  fre- 
quency as  required  by  the  adjustment  rules  of  Reference  6. 
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Figure  14.  T-38A  Pitch  Control  Dynamics;  Lag-Equalized  Pilot,  No  Bobweight; 
Effect  of  Pilot  Cain 


The  open-loop  0 -►  Fs  system  dynamics  with  the  fully  adapted,  lag  equalized 
pilot  are  shown  in  the  Bode  plot  of  Figure  14.;  the  stick-free  airplane 
dynamics  S/Fs(jio)  are  repeated  there  for  comparison.  Figure  14  also  shows  the 
closed  loop  system  amplitude  ratio  1 0/ 0 ( j to) | for  two  values  for  pilot  gain; 
the  value  Kp  = 156  lb/rad  was  the  one  suggested  in  Reference  4.  Recall  that  j 

the  pitch  attitude  command  0C  is  an  equivalent  command  due  to  turbulence  for 
conventional  VFR  flight;  it  could  also  be  a command  input  with  a flight 
director  system  or  an  "internally  generated"  command  by  the  pilot  himself. 

Actually,  0C  is  nearly  immaterial  to  the  PIO  analysis.  The  Important  criterion 
is  to  determine  the  nature  of  the  closed  loop  responses  i£  the  pilot  elects 
(or  is  forced)  to  execute  precision  control  of  0 — the  existence  or  non- 
existence of  command  or  disturbance  Inputs  is  irrelevant  to  this  purpose;  it 

is  an  analytical  convenience  to  "invent"  0 . 

c 


•>6 


The  closed  loop  response  ratio  |0/Oc(jw)|  shown  in  Figure  14  may  be 
used  to  compute  the  az  power  spectral  density  front  the  equation 
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Z(So  “ 4xM6e  Mwulz 
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I—  ut2 
T0  2 SP 


= (0. 5) 2 (ft/sec2)2  - sec 


Tt  is  assumed  that  the  pilot's  location  forward  of  the  e.g.  is  the  same  as 

the  effective  bobweight  location  (£  = 12.83  ft.  Reference  19).  In  view  of 

the  uncertainties  surrounding  0,  and  the  choice  of  an  appropriate  turbulence 

model,  assume  that  ^w^Wg^^  as  s'lown  In  Figure  15.  The  resulting  a7^ 

power  spectral  densities,  normalized  bv  no!  /10,  are  shown  in  Figure  lb  for 

g 

two  K values  with  the  lag-equalized  pilot  model.  These  two  spectra  are 
repeated  in  Figure  17  with  a linear  frequency  scale  for  ease  in  Interpreta- 
tion. Values  of  oa  are  tabulated  for  two  realistic  levels  of  cw  ; the 
index  v of  subjective  predictability  is  also  given  for  each  case.  The  values 
of  oaZp  were  obtained  by  direct,  approximate  integration  of  the  normalized 
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Figure  15.  Simplified  Turbulence  Model 
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Based  on  the  criterion  v>  r 0.3,  it  appears  that  az^(t)  is  sufficiently 
resonant  to  be  subjectively  predictable  with  both  pilot  gains — although  the 
case  for  = 98  is  marginal.  The  resonant  frequency  ' 3. 5-3. 6 rad/sec. 

The  third  case  shown  in  Figure  17  is  intended  to  simulate  the  situation 
where  the  pilot  drops  the  (typical)  lag  equalization  and  reverts  to  the  more 
primitive  proportional  control  of  0 with  a time  delay.  That  is 
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This  might  occur,  for  example,  when  the  pilot  becomes  "tense"  or  overly 
concerned  about  pitch  errors.  It  could  also  occur  when  the  pilot  does  not 
track  0,  at  all,  but  tracks  only  0;  there  is  some  basis  for  such  a model,  but 
it  is  beyond  the  scope  of  present  interest  and  will  not  be  considered  fur- 
ther. The  open  and  closed  loop  0 -*■  F^  dynamics  for  this  pilot  model  are 
shown  in  Figure  18.  Note  that  the  pitch  attitude  loop  closure  is  unstable 
for  Kp  = 156;  it  is  marginally  stable  (gain  margin  = 1 db)  for  = 98.  The 
power  spectrum  of  az(t),  for  the  rectangular  4>w  w (uj)  spectrum,  is  shown  in 
Figure  19  and,  as  previously  indicated,  in  linear  form  in  Figure  17.  It  is 
clear  from  Figure  17  that  the  dominant  effect  of  reducing  lag  equalization  is 
to  markedly  increase  the  resonance  and  resonance  frequency  of  uncontrolled 
a„  (t) . Figure  17  indicates  that  with  no  equalization  the  resonance  fre- 


quency u = 5.6  rad/sec. 
K 


Figure  18.  T-38A  Pitch  Control  Dynamics;  Non-Equalized  Pilot; 

No  Bobweight 
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Figure  19.  PSD  of  az^;  No  Bobweight;  Non-Equalized  Pilot 


It  is  clear  that  with  realistic  choices  for  Y^(ju))  the  pilot-felt 

acceleration  will  be  resonant  and  subjectively  predictable.  The  resonant 

frequency  w appears  to  be  very  sensitive  to  the  specific  form  chosen  for 
k 

Yp(jw);  it  is  reasonable  to  suspect  that  with  nominal  pilot  dynamics,  the 
resonant  frequency  will  be  much  smaller  than  that  predicted  for  off-nominal 
pilot  dynamics  resulting  from  the  reduction  or  elimination  of  lag  equaliza- 


The  implications  of  these  results  to  PIO  are  apparent  from  inspection  of 
— total  phase  of  the  az^  -*•  Fs  system — shown  in  Figure  20.  Two  cases  are 
shown  corresponding  to  no-  and  full-bobweight  (K  = 0 and  2,  respectively). 

D 

The  pilot's  time  delay  was  chosen  as  0.25  seconds  in  both  cases.  For  the 

small  acceleration  case  for  which  3 0,  it  can  be  seen  that  az^  -*■  F 

is  quite  stable  if  the  acceleration  loop  is  closed  such  that  it  resonates  at 

frequency  <*>  = 3.5  - 3.6;  this  is  the  case  for  the  fully-adapted,  lag  equalized 

K 

pilot  with  reasonable  gain.  The  phase  margin  is  about  80  degrees.  For  the 

non-equalized  pilot,  however,  the  phase  margin  is  markedly  decreased;  for  the 

case  of  K =98  and  w = 5.6,  it  can  be  seen  that  the  phase  margin  is  about 
P F 

20  degrees.  Thus,  it  is  plausible  that 
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Figure  20.  T-38A:  Total  Phase  Angle  of  the  Acceleration  Control 

Loop;  Effect  of  Bobweight 


1.  With  "normal"  pilot  dynamics,  and  for  initially  small  az^(t),  the 
control  of  9 can  lead  to  a7  resonance. 


2.  If  the  pilot  attempts  to  track  the  az^  cue,  the  resulting  az^ 

->■  F loop  dynamics  will  be  stable  and  az  amplifications  are  not 


likely  to  result. 


3.  With  a non-equalized  pilot,  the  control  of  9 can  lead  to  very 
resonant,  but  stable,  az^  responses. 
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4. 


- 


The  resonant  frequency  co  due  to  control  of  0 with  a non-equal i zed 

pilot  will  be  substantially  greater  than  that  resulting  from  pilot- 

adapted  lag  and  the  az  -*■  F loop  will  be  marginally  stable  at  the 

p s 

larger  wr. 

5.  The  probability  is  greatly  increased  that  az^(t)  will  become  of 
sufficient  amplitude  to  excite  the  bobweight  dynamics  when  the 
pilot  is  non-equalized . 

Of  course,  bobweight  dynamics  can  also  be  excited  with  the  fully-adapted 
pilot;  the  values  for  oa  shown  in  Figure  17  (for  Tj  = 1.2)  suggest  that 
it  would  be  a comparatively  rare  occurrence  when  az^  greater  than  1 g 
occurs . 


If,  for  whatever  reason,  the  amplitude  of  az^  becomes  sufficient  to 
excite  bobweight  dynamics,  then  further  tracking  of  0 will  result  in 

1.  A large  shift  in  w to  higher  frequencies  (typically  from  about  3.5 
to  about  8 rad/ sec). 

2.  Closed  loop  instabilitv  of  a.,  -*  F . 

^p  s 

This  appears  to  be  true  even  when  lag  equalization  is  used  provided  that  the 
pilot's  adaptation  is  with  respect  to  the  airplane  dynamics  corresponding  to 
Rj,  = 0.  This  situation  is  shown  on  the  Bode  diagrams  of  Figure  21;  open  and 
closed  loop  Bode  diagrams  are  shown  for  two  values  of  K with 


1 


Yp  (s) 


1 


as  before.  The  value  K = 15b  is  approximately  equal 

closed  loop  damping  rat io;  however,  even  with  K =98 

p 

highlv  resonant  with  .o  = 9.2  rad/sec.  Bv  comparison 

K 

it  is  clear  that  a.,  (t)  would  qualify  as  subjectively 


to  the  value  for  zero 
the  closed  loop  is 
with  previous  results 
predictable  for  these 


' 


Figure  21.  T-38A  Pitch  Control  Dynamics;  Lag- Equalized , Unadapted 

Pilot;  Full  Bobweight 


cases.  Inspection  of  <j>(  j w)  from  Figure  20  confirms  that  az  -*•  F would  be 

p s 

unstable  with  this  form  of  Y (s)  with  an  active  bobweight;  it  is  clear  that 

P 


closed  loop  instabilities  in  az  -*•  Fg  will  occur  for  greater  than  about 


8.5  rad/sec. 


Figure  20  bears  further  comment.  The  only  significant  difference 

between  the  az  /F  ( j 'o''  dynamics  with  and  without  the  bobweight  is  due  to  the 
P s 

short-period  mode.  The  phase  components  of  <t>(jw)  due  to  the  short-period 
with  and  without  the  bobweight  are  shown  in  Figure  20;  it  is  clear  that,  for 
a given  f requency  of  closure , az^  -*■  Fg  is  actually  more  stable  with  than 
without  the  bobweight  due  to  the  sizable  reduction  in  phase  lag  resulting 
from  the  bobweight  feedback  (chiefly,  the  stick-free,  short-period  damping 
ratio).  Thus,  by  present  theory,  if  the  pilot  attempts  to  track  a2  starting 
with  small  motions  and  no  effective  bobweight  and,  in  the  process,  creates 
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« 


a sufficient  to  excite  bobweight  dynamics,  the  effect  will  be  to  Increase 
the  az  -*•  F loop  damping.  If,  however,  the  pilot  reverts  to  the  0 * F 
loop  while  the  bobweight  is  still  excited,  he  effect  on  0 -*■  F loop  stability 


is  disastrous. 


The  above  analysis  is  offered  as  an  explanation  for  how  Type  1 Pit)  was 

initiated  with  the  T-38A.  The  only  significant  nonlinearity  involved  is 

control  system  friction;  its  effect  is  to  mask  bobweight  effects  on  stick- 

free  dynamics  below  a threshold  level  of  normal  acceleration.  Once  the 

pilot  begins  to  track  large- amplitude  a..  , however,  the  breakout  force  will 

P 

no  longer  be  of  primary  significance. 


There  are  other  control  system  nonlinearities  that  are  important  to 
consider  when  the  oscillations  become  of  large  amplitude.  Chief  among  these 
is  probably  the  combined  effect  of  nonlinear  control  gearing  and  control 
system  hysteresis  (Reference  19).  The  sinusoidal  describing  function  for 
this  is  given  in  Reference  19  and  is  reproduced  in  Table  3 for  the  case  of 


oscillations  about  trim.  The  a„  control  system  may  be  modeled  as  indicated 


in  Figure  22a.  It  can  be  shown  that  an  equivalent  form  for  the  indicated 
system  dynamics  is  that  depicted  in  Figure  22b;  this  is  the  form  required 
for  the  determination  of  system  limit  cycles  via  the  gain-phase  analysis. 
By  inspection. 


TABLE  3.  SINUSOIDAL  DESCRIBING  FUNCTION  FOR 
T-38A  CONTROL  SYSTEM 


| Fs  | 1 

b max 

(lb) 

1 

Gain  |n| 

Phase  ^ N 

(deg) 

(deg/lb) 

(db) 

5 

.10 

-20 

-32 

10 

.18 

-14.9 

-12 

20 

.32 

-9.9 

-13 
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Figure  22.  T-38A:  A Model  for  Large-Amplitude  Acceleration 

Control  Dynamics 


L(jw)  = G(j  w)  [KD+P  (j  oi)  ] = dynamics  of  the  linear  portion  of  the  system. 

D 

N = feel  system  describing  function  (Table  3). 

G(jw)  forward  path  dynamics  of  the  linear  control  system  and  airframe; 

G(j“)  = [TwVTfoiTlfs^isTTW  g/lb  ; (the  static  sain  ls  0,508  g/lb)‘ 

Kg  = bobweight  gain  for  large  amplitude  motion;  for  the  original  T-38A 
control  system  Kg  = 2 (lb/g). 

. „ -t  iu  = pilot  dynamics  for  a7  control;  it  is  assumed  that 

P(ju))  = K e aJ  J zp 

Si 

x = 0.25  sec. 
a 


The  gain-phase  plots  for  both  L(jco)  and  -1/N  are  shown  in  Figure  23; 
five  values  of  the  pilot  gain  K^  were  used  for  the  computation  of  L(jw). 

It  can  be  seen  that  a stable  limit  cycle  can  exist  with  amplitude  and 
frequency  dependent  upon  pilot  gain;  typical  results  are  shown  in  Table  4. 
The  acceleration  amplitudes  shown  in  Table  4 were  computed  directly  from  the 


a,  / F (iui)  transfer  function. 

An  O J 


If  aerodynamic  nonlinearities  at  high  angles 


of  attack  were  included  in  this  analysis,  then  the  unrealistic  accelerations 
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Now,  according  to  the  P10  theory  described  previously  in  tills  report,  a 
Type  1 P10  can  develop  provided  the  accelerations  resulting  from  pitch 
attitude  control  are  sufficiently  periodic  (and  it  was  shown  above  that  they 
are,  in  general)  and  that  the  resonance  frequency  is  approximately  equal 
to  or  greater  than  w,  . It  was  shown  above  that  w is  highly  variable 
depending  upon  pilot  adaptation  and  whether  or  not  the  bobweight  dynamics 
are  excited.  For  those  conditions  believed  to  be  related  to  the  pre-PIO 
phase,  however,  the  l'  -*•  F loop  analyses  suggests  that  will  lie,  typi- 
cally, between  5.6  and  9.5  rad/sec.  We  conclude,  therefore,  that  it  large 
amplitude  motions  are  initiated  and  If  the  pilot  attempts  to  track  az  , 
then  a stable  limit  cycle  will  be  very  likely.  The  control  system  hysteresis 
at  large  amplitudes  is  sufficient  to  sustain  the  limit  cycle. 


The  third  and  final  necessary  condition  for  Type  I P10  was  that 


aZp 

6 


(jWR> 


> 0.012 


g/degree /sec 


This  ratio  for  the  T-38A  is  given  in  Table  5 for  various  possible 
Clearly,  the  amplitude  criterion  is  satisfied. 


TABLE  5. 
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(rad /sec) 
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(g/deg/sec) 
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.196 
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It  is  concluded  that  Type  I PIO  is  possible  for  this  airplane,  that  the 
incipient-PIO  frequency  will  be  variable  (from  about  5.6  to  9.5  rad/sec), 
and  that  the  fully-developed  PIO  frequency  will  be  that  of  a stable  limit 
cycle  occurring  in  the  az^  -*■  Fg  loop  of  frequency  wpjQ  = 6.1  to  7.2  rad/sec. 
Present  theory  substantiates  the  basic  conclusions  of  previous  investiga- 
tions of  the  system-oriented  causes  of  the  problem  (notably  References  1,  4, 
and  19);  the  present  explanation  of  the  phenomenon  is,  however,  qualitatively 
different. 


Type  II  PIO  Assessment — T-38A 


From  the  T-38A  data  summarized  in  Table  2,  it  can  be  seen  that  for 

large  amplitude  control  or  disturbance  inputs  the  stick-free  short-period 

damping  ratio  C = 0.1;  by  the  simplified  criterion  for  subjective  pre- 
sp 

dictability  of  az^,  conclude  that  the  pilot  can  recognize  the  periodicity  of 

the  dominant  mode  response  for  abrupt  inputs  and  he  may  attempt  to  close  the 

a_  -*■  F loop  at  = 9.8  rad/sec.  This,  however,  can  result  in  a large- 
Ap  S R 

amplitude  limit  cycle  as  shown  above.  It  can  be  seen  from  Table  5 that  the 
amplitude  criterion  is  satisfied  at  = 9.8  rad/sec.  It  is  therefore 
apparent  by  present  theory  that  the  T-38A  with  the  original  control  system 
was  susceptible  to  Type  II  PIO. 


PIO  Assessment — Modified  Control  System 

A number  of  modifications  were  made  to  the  original  T-38A  control 
system  in  an  attempt  to  cure  the  airplane's  PIO  problems.  The  most  signi- 
ficant of  these,  so  far  as  present  theory  is  concerned,  were 

1.  Increased  feel  spring  rate  at  small  stick  deflections. 

2.  Decreased  bobweight  gain  from  about  2.0  to  1.0  lb/g  measured  at 
the  control  stick  grip. 

The  intent  of  these  modifications  was  to  reduce  the  pilot's  tendency  toward 
"overcontrol"  and  to  reduce  the  large  variation  in  stick-free  dynamics  from 
small  to  large  levels  of  normal  acceleration,  while  maintaining  approxi- 
mately the  same  level  of  stick  force  per  g (References  1 and  4). 
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It  is  clear  from  Table  2 that  with  the  modified  control  system  there  is 
no  longer  a substantive  difference  between  the  stick-fixed  and  the  stick- 
free  dynamics.  Consequently,  the  6 -*■  Fg  analyses  previously  dis-'ussed  for 
the  case  of  small  amplitude  motion  (Kg  = 0)  are  also  approximately  valid  for 
motion  of  any  amplitude  with  the  modified  control  system.  It  can  be  con- 
cluded, then,  that  while  closed  loop  resonance  can  still  occur,  the  chances 
for  it  are  greatly  reduced.  The  most  important  consideration,  however,  is 
that  for  realistic  pilot  dynamics  (with  or  without  fully-adapted  equaliza- 
tion) the  resonant  frequency  will  be  less  than  that  required  for  instability 
of  the  acceleration  loop  az^  -*■  Fg  when  control  system  nonlinearities  are 
ignored.  In  such  cases.  Type  I PIO  is  impossible  by  present  theory. 


The  modifications  made  to  the  control  system  do  not  appear  to  have  had 
a substantial  effect  on  the  large-amplitude  control  system  nonlinearities. 
Figure  24  is  a gain-phase  plot  of  the  linear  and  nonlinear  portions  of  the 
a^  Fg  system  dynamics  for  the  modified  system.  The  only  difference 
assumed  between  this  analysis  and  the  one  summarized  in  Figure  23  is  that 
the  bobweight  gain  is  1.0  lb/g  rather  than  2.0.  It  is  clear  from  the  gain- 
phase  portrait  that  if  large  amplitude  resonant  motions  can  be  initiated  for 
any  reason  with  to  greater  than  about  6.0  rad/sec,  and  if  the  pilot  attempts 
to  track  az^,  then  a stable  limit  cycle  can  still  occur  with  the  modified 
flight  control  system. 


A recent  PIO  occurrence  with  the  modified  T-38A  does  appear  to  corrobo- 
rate the  present  theoretical  conclusion.  This  incident  occurred  with  Air 
Training  Command  Aircraft  SN  64-3253  in  November  1975.  The  flight  condition 
was  M = 0.87  at  an  altitude  of  15,000  feet.  The  aircraft  was  successfully 
recovered.  Other,  similar,  problems  with  the  T-38A  -have  apparently  occurred; 
these  were  the  subject  of  Reference  21.  Subsequent  investigation  determined 
that  a f lap-stabilator  interconnect  cable  had  broken  prior  to  the  PIO  and 
that  the  probable  effect  of  this  breakage  was  to  increase  the  pitch  control 
gain  by  a factor  of  as  much  as  3.0  from  the  normal  value.  The  resulting 
pitch  sensitivity  led  to  the  excursions  in  acceleration.  The  instructor 
pilot  was  at  the  controls  when  the  PIO  occurred.  If  it  is  assumed  that  the 
pilot  was  fully  adapted  with  respect  to  the  normal  aircraft  dynamics,  then 
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the  previous  analysis  with  the  lag-equalized  pilot  (1/Tf  = 3.2),  shown  in 
Figure  14,  is  approximately  equivalent  to  this  situation  as  a model  for 
pitch  loop  dynamics.  (Note  that  the  no-bobweight  case  with  the  unmodified 
control  system  is  dynamically  similar  to  the  full-bobweight  case  with  the 
modified  system.)  If  a realistic  value  for  is  98  lb/rad,  then  it  can  be 
seen  from  Figure  14  that  a factor  of  3 increase  in  open  loop  gain  (i.e.,  9.5 
db)  will  yield  an  unstable  pitch  loop  closure;  the  resulting  motions  in  0 
will  be  a divergent  oscillation  with  frequency  approximately  equal  to  7.4 
rad/sec.  Inspection  of  <J> ( j to ) (Figure  20)  indicates  that  the  acceleration 
loop  az  -*■  F will  also  be  unstable  when  closed  at  this  frequency.  Thus, 
linear  analysis  indicates  that  the  non-gain-adapted  pilot  will  create  large 
amplitude  motions  due  to  the  control  system  linkage  failure.  Figure  24 
suggests  that  a large  amplitude  limit  cycle  would  be  the  probable  steady 
state  form  of  the  motion.  It  was  Indicated  that  in  this  PIO  encounter  the 
pitch  oscillations  continued  for  several  seconds  and  that  the  post-flight 
accelerometer  readings  were  plus  5.2  and  minus  6.0  g's.  by  comparison  with 
past  T-38A  PIO  case  histories,  this  suggests  that  a stable,  large-amplitude 
limit  cycle,  resulting  from  az^  loop  closure,  is  a reasonable  explanation 
for  tlie  facts  surrounding  this  incident.  Note  the  consistency  between  the 
g's  obtained  in  the  PIO  and  those  indicated  in  Figure  24  or  Table  4 for  a 
factor  of  three  increase  in  nominal  pilot  gain;  i.e.,  the  g level  decreases 
with  increasing  gain. 


Type  l l PIO  would,  by  present  theory,  appear  to  have  been  eliminated 
since  the  dominant,  stick-free  mode  damping  ratio  4 ' = 0.28  and  is  greater 

than  the  criterion  value  of  0.2.  It  is  therefore  concluded  that  if  PIO  is 
to  be  initiated  with  the  modified  control  system,  then  it  must  result  from  a 
m isadaptat ion  by  the  pilot  (excessive  gain  or  lag  drop-out  in  the  control  of 
pitch)  together  with  a concern  for  load  factor  control. 
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Figure  24.  T-38A:  Gain-Phase  Diagrams  for  the  Nonlinear 

aZp  * Fs  System  Dynamics;  Modified  Control  System 

Comment  on  Pilot  Time  Delay 

It  was  indicated  in  Section  III  that  with  a physiologically  acceptable 
model  for  pilot  dynamics,  it  is  not  necessary  to  resort  to  the  assumption 
that,  in  P10,  the  pilot  is  synchronous  (i.e.,  has  no  time  delay  or  equaliza- 
tion). It  has  been  assumed  that  a small  pilot  delay  exists  in  the  az  ■+ 
loop  and  that  a reasonable  value  for  this  delay  might  be  0.25  seconds  for 
want  of  better  data. 

In  fact,  present  theory  requires  a minimum  value  of  T in  order  to 
support  a P 1 0 . This  can  be  seen  with  the  present  example.  In  Figure  24  the 
dashed  curve  is  the  gain-phase  plot  for  L(ju>)  with  K = 8 lb/g  and  T = 0.17. 

<1  u 

The  two  curves  (L  and  — 1 /N)  are  tangent  at  8 s u>  s 0 but  are  otherwise  non- 
intersecting. The  conclusion,  therefore,  is  that  to  obtain  PIO  by  present 
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theory  T must  be  greater  than  0.17  seconds.  The  value  T = 
a a 

selected  to  give  uniformly  reasonable  PIO  amplitude  and  frequencies  for  the 
available  case  history  data.  This  point  will  be  noted  again  in  the  YF-12 
analysis . 

C.  THE  YF-12  (M  = 0.77,  ALTITUDE  = 25,000;  REFUELING  CONDITION) 

Reference  11  documents  and  analyzes  PIO  of  two  distinct  types  that  are 
suspected  to  have  occurred  with  this  aircraft,  SAS-on,  in  the  same  flight 
c ond it  ion. 

The  more  serious  of  these  PIO  was  indicated  to  have  occurred  at  least 
twice.  It  is  characterized  by  oscillations  in  normal  acceleration  of  about 
' 2 g at  a frequency  of  about  1/2  cps.  Reference  11  implicates  saturation  of 
the  pitch  damper  due  to  large  amplitude  motions  as  being  central  to  its 
development . 

The  second  PIO  mode  was  indicated  to  be  of  about  1 cps  in  frequency 
with  an  amplitude  of  about  ±1/4  g.  Reference  11  states  that  this  oscillation 
is  rather  commonly  encountered  during  refueling  and  that,  while  it  is 
annoying  to  the  pilot  and  quite  beyond  his  ability  to  control,  it  is  other- 
wise benign.  Its  origin  was  attributed  to  coupling  between  structural 
dynamics  and  rigid  body,  short-period  dynamics. 

Reference  11  presents  an  analysis  of  both  these  PIO  modes  based  on 
dynamic  properties  of  the  closed  loop  pitch  attitude  control  system.  The 
pilot  model  used  is  a pure  gain,  following  Reference  19  and  others.  The 
results  of  Reference  11  are  essentially  qualitative;  they  serve  to  illustrate 
the  importance  of  SAS  amplitude  or  rate  limits  to  the  initiation  of  large- 
amplitude  oscillations.  The  authors  note,  in  addition,  that  their  analysis 
and  similar  past  analyses  of  the  PIO  problem  are  not  capable  of  PIO  pre- 
diction but  serve  instead  to  identify  causal  factors  after  PIO  is  known  to 
occur. 
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report  to  explain  both  PIO  inodes  of  the  YF-12,  thereby  confirming  both  the 
theory  and  the  principal  results  of  Reference  11.  This  is  only  possible  in 
part,  unfortunately,  due  to  analytical  complications  posed  by  the  SAS  non- 
linearities  and  to  the  almost  complete  lack  of  an  analytical  model  for  the 
effects  of  structural  dynamics  on  the  az  /F  transfer  function.  It  will  be 
shown  in  the  following  that  the  large  amplitude,  potentially  catastrophic 
PIO  is  rather  easily  predicted  by  present  theory.  The  "nuisance"  PIO  cannot 
be  quantitatively  confirmed  by  direct  application  of  the  theory  due  to  the 
lack  of  a suitable  model  for  structural  dynamics.  It  will  be  shown,  however, 
that  its  character  can  be  assessed  with  a "plausible"  model  for  structural 
effects  and  that  the  resulting  implications  to  pitch  attitude  and  pilot's 
control  response  dynamics  are  consistent  with  the  flight  test  results  pub- 
lished in  Reference  11.  In  fact,  the  analysis  offered  here  for  the  "nui- 
» sance"  PIO  appears  to  be  the  simplest  and  most  straightforward  method  by 

which  a theoretical  understanding  of  the  control  deflection  power  spectral 
density,  obtained  in  flight  tests  of  the  YF-12,  can  be  derived. 


The  dynamics  of  the  YF-12  control  system  are  complicated  by  a flexible 
fuselage  structure  and  by  SAS  nonlinearities.  It  is  assumed  that  an  appro- 
priate block  diagram  for  the  pitch  attitude  closed  loop  is  as  shown  in 
Figure  25.  The  rigid-body  modal  response  9 is  considered  to  be  the  appro- 
priate  SAS  input,  since  the  rate  gyro  is  placed  near  a point  of  zero  slope 
on  the  first-mode  fuselage  bending  curve.  The  pilot,  however,  is  assumed  to 
sense  the  combined  attitude  0 due  to  rigid-body  motions  and  first  mode 
bending. 


In  general,  the  analysis  of  such  a system  is  impossible  except  by 
computer.  Fortunately,  it  is  possible  to  ascertain  the  qualitative  nature 
of  the  pilot's  control  problem  by  examining  only  the  limiting  cases  of  very 
small  motion  amplitudes  for  which  the  SAS  is  completely  linear  and  large 
motions  for  which  the  SAS  may  be  considered  to  be  always  saturated  and. 


therefore,  nonfunctional.  Thus,  the  models  required  are  0/6^ 
in  the  SAS-on  and  SAS-off  conditions. 


and  az  /S 
P e 


First  Mode 
Bending 


Figure  25.  YF-12  Pitch  Attitude  Control  Model 


It  should  be  noted  that  the  YF-12  elevon  controls  are  much  more  compli- 
cated than  indicated  by  Figure  25.  The  differences  are  believed  to  be 
unimportant  for  an  analysis  of  manual  control.  At  most,  higher  frequency 
dynamics  associated  with  the  actuators  contribute  a small,  equivalent  time 
delay  of  about  0.06  seconds  to  the  0 loop  dynamics.  Actuator  dynamics  will 
be  neglected  in  the  remainder  of  this  analysis. 

The  first-mode  structural  bending  curve  is  reproduced  below  (Figure  26) 
from  Reference  11.  In  private  communications  with  the  authors  of  Reference 
11,  it  was  disclosed  that  this  curve  was  analytically  derived. 
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(From  Reference  11) 
igure  26.  YF-12  First  Mode  Fuselage  bending  Curve 
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The  stability  derivatives  given  in  Reference  11  are  related  to  those 
used  in  this  report  (all  of  which  are  consistent  with  the  notation  of 
Reference  20)  according  to  Table  6. 


TABLE  6.  STABILITY  DERIVATIVES  FOR  THE  YF-12 


This  Report 

(Reference  2) 
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UD  = 786  ft/sec 

£x  = 50  ft  (from  Fig.  26) 

M = 0.8 

h = 25,000  ft 

Airframe  Dynamics  (No  SAS;  Rigid  Body) 
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Airf rame  Dynamics  (With  SAS;  Ki^id  Body) 
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Airframe  Dynamics  (No  SAS;  First  Bending  Mode  Included) 
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From  Reference  20, 
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If  the  pilot's  displacement,  due  to  fuselage  bending,  from  the  line  of  zero 
bending  (Figure  26)  is  Zp,  then 

zp  = “ ^o^B 

zp  = - Mb 

where  the  pilot's  location  forward  of  the  node  of  the  first  bending  mode  curve 
is  S.G  = 19.3  feet.  If  aZp  is  assumed  to  result  only  from  the  rigid  body 
acceleration  plus  the  bending-induced  acceleration,  then 


262 (s2+2 ( . 114) (5. 35)s+5. 352  ] [s2+2 ( . 044) (11. 91)s+ll. 912 1 
[s2+2 ( . 376) (2. 01) s+2 . 012 ) [ s2+2 ( . 050) (15 . 7)s+15 . 72 ] 


f t/sec2/rad 


Unfortunately,  this  result  neglects  the  heave-induced  component  of  aZp  for  which 
there  is  insufficient  data  available  for  a reliable  estimation.  Furthermore, 
there  is  no  way  to  assess  the  error  that  neglecting  this  mode  might  entail.  For 
this  reason,  the  results  of  this  section  are  qualitative. 

Airframe  Dynamics  (With  SAS  and  First  Bending  Mode) 


-11 . 23 (s+. 79) (s+3.89) [s2+2(,10) (11 . 78) s+11 . 78? ] 

s (s+1 .478) [s?+2 ( . 692) (4 . 56) s+4 . 562 ][ s24 2 ( . 050) (15. 7)s+15. 72 ] 


rad/rad 


With  the  previous  assumption  that  bending-induced  heave  is  neglected,  obtain 
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Feel  System  Dynamics 


The  feel  system  dynamics  are  indicated  by  Reference  20  to  be  very 
nonlinear  as  follows: 


o There  is  a large  breakout  force  (about  5 pounds). 

o The  input  of  sinusoidal  stick  force  Fs(t)  produces  a sizable 
hysteresis  loop  in  the  control  stick  deflection  response. 


o 


The  control  gearing  ratio  ^ec/^e  nonlinear* 


The  series  combination  of  feel  system  hysteresis  and  nonlinear  control 
gearing  was  used  in  Reference  11  to  develop  a sinusoidal  describing  func- 
tion. This  is  reproduced  in  Figure  27. 


Figure  27. 


Describing  Function  for  YF-12  Control  System 


Type  I PIO  Assessment — YF-12 

The  assessment  of  the  YF-12' s potential  for  Type  I PIO  is  analogous  to 
that  for  the  T-38A.  The  pitch  dynamics  of  the  YF-12  qualitatively  change 
when  pitch  rate  becomes  sufficient  to  saturate  the  position-  and  rate- 
limited  SAS;  the  bobweight  and  control  friction  produced  similar  behavior  in 
the  case  of  the  T-38A. 

Figure  28  is  a Bode  plot  of  the  linear  B/Sg^  dynamics  for  the  YF-12 
with  and  without  the  SAS  and  including  flexibility  effects.  Also  shown  is 
the  computer-predicted,  nonlinear  frequency  response  taken  from  Reference  11 
for  the  case  where  the  SAS  is  assumed  to  be  position  limited  to  2.5  degrees 
authority  and  with  a rate  limit  of  12.6  degrees/sec.  It  is  clear  that  when 
the  SAS  limits  are  encountered  the  resulting  pitch  dynamics  are  given  to  a 
good  approximation  by  the  SAS-off  dynamics — a very  fortuitous  result!  It 
was  stated  in  Reference  11  that  variations  in  the  position  limits  o'.  SAS 
authority  do  not  greatly  affect  the  nonlinear  frequency  response;  the  rate 
limit  is  the  essential  SAS  parameter. 

Inspection  of  Figure  28  suggests  that,  to  a good  approximation,  0/ <5e^  is 
of  the  form  K/s  in  the  vicinity  of  the  probable  crossover  frequency  (3-4 
rad/sec)  when  the  SAS  is  on  and  not  limiting  (the  dashed  curve).  Thus,  an 
appropriate  pilot  model  for  the  control  of  pitch,  with  small  amplitudes  of  ti, 
would  be 

Y (s)  = K e * 

P P 

i.e.,  a gain  with  time  delay.  It  can  be  seen  from  Figure  28  that  the 
maximum,  permissible  gain  for  stability  of  the  pitch  loop  is  reduced  by  a 
factor  of  about  3.75  (11.5  db)  when  the  motions  become  of  size  sufficient  to 
saturate  the  SAS.  Also  shown  in  Figure  28  is  a closed-loop  pitch  response 
transfer  function  for  the  case  where  the  SAS  is  unsaturated  and  the  pilot 
gain  is  large.  It  is  clear  that  the  pitch  dynamics  can  be  resonant  with  a 
frequency  u>  s?  4.5  - 4.9  rad/sec;  this  is  also  the  most  likely  range  of 

maximum  values  for  u>  . 

K 
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Figure  28.  YF-12  Pitch  Attitude  Dynamics 


In  Figure  29  the  open-loop  transfer  function  is  plotted  for  the 
F loop  with  a linear  (i.e.,  unsaturated)  SAS; 


Fs  Dynamics;  SAS-On 


Figure  29.  YF-12  aK  " 

v , . „ -.25s  azP  , * 

Yql  (s)  - Kp  e (») 

The  total  open- loop  phase  angle  Is  shown  for  two  cases:  one  will)  no  phase 
lag  due  to  control  and  feel  system  dynamics  and  one  with  a time  delay  of 
0.16  seconds  (0.06  seconds  due  to  the  actuator  lags  and  0.10  seconds  due  to 
the  feel  system  lag).  In  both  cases  the  pilot  delay  is  assumed  to  ho  0.25 
seconds . 


it  Is  apparent  from  Figure  29  that  az  ^ *■  F system  Instability  is 
possible  at  a resonant  frequency  w ^ *s  4.5  - 4.9  rad/sec  when  the  loop 
del. iv  Is  about  0.  i2  seconds.  This  could  occur  when  the  control  system 
delays  are  approximately  equal  to  or  greater  than  0.07  seconds — not  an 
unreasonable  value  for  the  YF-12  control  system. 


t i me 


To  assess  Che  large-amplitude  character  of  terminal  oscillations 
resulting  from  any  Instability  of  the  linear  system  for  small-amplitude 
motions,  it  is  necessary  to  examine  the  gain-phase  plots  of  the  serialized 


linear  and  negative  reciprocal,  nonlinear  portions  of  the  az ^ ♦ F^  system, 


as  was  done  for  the  T-38A.  For  the  YF-12,  we  assume  that  with  large- 


amplitude  motions  the  SAS  effects  on  az  /dp^  dynamics  may  be  neglected  and 


the  appropriate  transfer  function  Is  that  for  no  SAS.  The  linear  part  of 
the  az  -►  F loop,  based  on  previous  assumptions,  is  then 


. , . ..  -.25s  zp  . . 

l.(s)  - k e — (s) 


The  bode  plot  for  l.(Jio)  Is  shown  in  Figure  30.  The  nonlinear  part  MvA)  is 
the  describing  function  for  feel  and  control  system  nonlinear  1 1 ies  shown  In 
Figure  27.  The  gain-phase  plots  of  l.(|u)  and  -1/N  are  shown  in  Figure  11 


for  a pilot  gain  of  - 16.6  lb/g.  It  is  clear  that  a stable  limit  cycle 


can  be  obtained  at  a frequency  of  approximately  3.1  rad/sec;  the  control 
amplitude  is  dependent  on  K^;  the  value  - 16.6  lb/g  gives  |<S0^|  ■>  8.0 


degrees . 


Note,  also,  that  Figure  31  indicates  that  L(|io)  and  -1/N  are  tangent 


but  otherwise  non- intersect ing  for  K «*  16.6  when  t “ 0.22  (the  dashed 


•urve).  This  result  is  not  too  sensitive  to  K . This  is  turther  evidence 

a 


to  support  the  assumption  that  r " 0.25  is  a generally  valid  delay  for 

Jl 


study  ol  t lie  normal  acceleration  dynamics. 
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Thus,  the  estimated  peak-to-peak  oscillations  In  the  terminal  limit  cycle 
are  A g.  These  predictions  of  initial  frequency,  terminal  frequency,  and 
I’U'  amplitude  closely  match  the  YF-12  large-amplitude  l’lO  time  history 
published  in  Reference  20  and  reproduced  here  as  Figure  12. 


(From  Reference  l 1 ) 

Figure  12.  YF-12:  Large-Amplitude  1’ 10  Time  History 


The  large  amplitude  1M0  can  also  he  Initiated  hv  large  amplitude,  high 
gain  tracking.  It,  due  to  turbulence  or  large  amplitude  control,  the 
motions  become  large  enough  to  saturate  the  SAS  then  it  is  clear  from 
Figure  28  that  the  pitch  loop  can  be  destabilized  unless  t Ho  pilot  reduces 
tils  0- loop  gain  by  a factor  of  about  4.  It  the  SAS  is  saturated  then  it  can 


ho  seen  from  Figure  10 
t requencv  greater  than 


that  the  a_  * F 
Zp  s 


about  1 rad /sec. 


loop  is  unstable  when  tuned  t o anv 
It  Is  therefore  apparent  that  t ho 


large  amplitude  YF-12  Pit'  can  result  from  either  a^  * F loop  Instability 

when  the  0 ► F loop  is  driven  to  resonance  with  unsaturated  SAS,  or  from 
s 

a » F loop  instability  due  to  SAS  saturation  with  mlsadapted  pilot  gain, 
s 

In  either  case  the  terminal  state  ot  motion  is  a stable  limit  cvcle. 


The  small  amplitude  YF-12  P10  is  more  difticult  to  analyze  in  a del  ini- 
rive  manner.  it  is  not  cleat  from  a priori  consider at  ions  that  it  is  a t t tie 
1*10  as  defined  In  t'hapter  11. 


8 a 


Figures  33  and  34  are  adapted  from  Reference  11  to  show  the  flight  test 
measured  power  spectral  densities  of  0 and  |5e}.  during  refueling  while  the 
YF-12  is  actually  connected  with  the  tanker.  The  implications  of  these 
data  to  human  pilot  dynamics  offer  additional  understanding  of  the  PIO 
phenomenon . 


If  it  is  assumed  that,  during  those  periods  when  the  small  amplitude 
PIO  exists,  the  SAS  effects  on  pitch  dynamics  are  negligible  (e.g.,  if  it 
is  saturated),  then  from  the  Bode  plot  of  Figure  28  it  can  be  seen  that 


0/6e  is  approximately  equal  to  a pure  gain  plus  delay  in  the  region  of 
2 rad/sec.  Thus,  for  precision  tracking,  an  appropriate  model  for  pilot 
dynamics  might  include  a sizable  lag  equalization.  Open  and  closed  loop 
Bode  diagrams  are  shown  in  Figure  35  for  the  pilot  model 


I . 

(.  From  Re f e rence  1 1 ) 

Figure  J3.  YF-12:  PSD  of  Pitch  Attitude  During  Refueling 


i 


Ihe  correspond  Inn  power  spectrum  tor  0 will  obviously  be  qualitatively  con- 
sistent with  the  t light  test  result  shown  in  F inure  13;  the  calculation  for 
this  will  not  be  shown,  however.  It  is  noteworthy  that  the  replication  of 
these  data  by  analysis  appears  to  require  substantial  pilot  lan  equalization; 
otherwise,  the  low  frequency  peak  in  the  0 PSP  cannot  be  generated.  This 
result  is  independent  ol  whether  or  not  the  SAS  is  saturated;  it  is  suf- 
ficient that  the  pilot  model  contain  the  lag.  The  lag  is  easier  to  justify 
on  empirical  grounds  (viz.,  Reference  b)  when  the  SAS  is  saturated.  The 
correspond ing  PSP  of  the  pilot's  stick  force  deflection  Se  , due  to  pitch 
attitude  control,  is  shown  in  Figure  34;  it  is  apparent  that  the  pitch  loop 
closure  accounts  tor  the  low  frequency  peak  in  It  does  not  appear 

that  the  3\SiS  peak  at  l cps  can  be  explained  on  the  basis  of  pitch  control. 

It  has  been  noted  that  the  estimation  ot  a.  /F  dynamics  is  not  really 

P s 

possible  with  the  information  provided  by  Reference  11.  In  Figure  10,  a 

sketch  ot  "plausible"  a-,  /F  dynamics  is  shown  (the  dashed  curves)  which  are 

p s 

realistic  in  view  ot  known  or  suspected  structural  dynamic  effects.  The 

ground  rules  for  selection  of  these  "plausible"  dynamics  were  that 

o They  produce  substantial  structural  coupling  at  1 cps 

rite  a.  *■  F closed  loop  system  must  be  resonant  at  about  l cps  when 
P » 

the  pilot’s  gain  is  chosen  to  make  the  acceleration  loop  crossover 
frequency  approximately  equal  to  the  resonant  frequency  of  the  pitch 
i oop 

o 1'hev  are  consistent  with  the  large-amplitude  prediction  of  I'lO 

rite  closed  loop  a-  * F system  dynamics  are  shown  in  Figure  30  for  the  case 
p s 

where  the  pilot's  gain  is  selected  to  make  the  crossover  frequency  equal  to 
the  pitch  loop  resonance  frequency  (1.5  rad/sec).  The  corresponding  control 
del  lection  power  spectrum  was  computed  and  is  shown  in  Figure  34. 

It  we  now  assume  that  the  total  control  PSD  is  equal  to  the  sum  of  the 
components  due  to  & and  az  control  then  the  result  is  shown  by  the  dashed 
curve  in  Figure  34.  tlearly,  this  result  is  qualitatively  and  quantita- 
tively consistent  with  the  flight  test  result.  There  is  no  other  method  to 
explain  these  results  that  is  obvious  to  this  writer.  There  are  two  note- 
worthy features  ot  these  postulated  0 and  a2  loop  dynamics.  First,  at  the 


4 


frequency  of  resonance  of  the  az  loop  (1  cps)  the  pitch  attitude  response 

is  negligible;  thus,  there  is  no  loop  coupling  from  a,  » F to  u)  * F . 

p s s 

Second,  if  the  pilot's  acceleration  loop  gain  were  selected  to  make  the 
crossover  frequency  equal  the  az  resonance  frequency  (1  cps)  then  a closed 
loop  instability  could  result;  this  apparently  does  not  happen.  We  haven't 
the  data  to  determine  why  it  does  not.  One  possibility,  of  course,  is  that 


the  az  /F 
P s 


closure  at 


dynamics  (if  these 
1 cps.  This  seems 


were  known  in  detail)  would  support  a stable 
Implausible.  Another  possibility  is  that  the 


amplitude  of  az  resonance  at  1 cps  is  simply  too  small  to  permit  the  pilot 
to  subjectively  predict  az^  or  is  too  small  to  be  of  overriding  concern  to  him. 
rills  possibility  cannot  be  evaluated  without  a reliable  model  for  structural 


dynamics . 


it  is  concluded  that  by  present  theory  the  1 cps  oscillation  of  the  YF-12 
documented  in  Reference  11  is,  indeed,  a PIO  by  the  definition  of  this 
report.  It  is  a lightly-damped  resonance  in  acceleration  that  results  from 
attempts  by  the  pilot  to  stablize  the  normal  acceleration  motions  resulting 
from  closed  loop  dynamics  of  the  pitch  loop. 


Type  II  IMP  Assessment --YF - 1 2 

According  to  the  simplified  theory.  Type  II  F10  is  not  possible  in  the 

refueling  configuration.  The  stick-free,  dominant  mode  damping  ratio  is 

7.  ■ 0. 37fr  (SAS-off);  this  is  greater  than  the  criterion  value  of  0.2. 

sp 

b.  CALSPAN  P10  IN-FLir.HT  SIMULATIONS 

Reference  17  documents  an  in-flight  simulation  with  the  variable 
stability  NT-33A  of  approximately  150  combinations  of  short-period  frequency, 
damping  ratio,  and  speed.  Stick  force  per  g and  elevator-to-st ick  deflec- 
tion gearing  were  pilot  selected.  The  purpose  of  these  experiments  was,  in 
part,  to  examine  the  suitability  of  the  parameter 

2i  u> 

SP  sp 

L 


as  a measure  of  Pit)  susceptibility.  It  had  been  proposed  that  values  of 
this  ratio  less  than  one  were  a prerequisite  for  PIO,  and  that  values 
greater  than  one  would  be  sufficient  to  preclude  the  development  of  P10. 

This  criterion  originated  from  the  notion  that  in  P10  the  pilot  is  a syn- 
chronous controller  of  pitch  attitude  (discussed  previously  in  Section  11 1 . 

Sustained  P10  was  obtained  for  only  six  of  the  short-period  dynamic 
cout igurat ions  tested.  One  of  those  ft  light  58  0 violated  the  above  cri- 
terion. Sustained  PIO  was  defined,  in  Reference  I',  to  correspond  to  a P10 
rating  (P10R>  ot  -»  or  greater  according  to  the  I’lOK  scale  devised  (apparently! 
at  McDonnell  Aircraft  during  t light  tests  ot  the  K-s.  1'his  scale  is  shown 
in  Figure  In.  These  six  coni'  igurat  ions  are  summarized  in  Table  7. 

I'he  theory  of  this  report  may  be  easily  applied  to  these  cent' igurat  ions, 
"his  analysis  will  not  he  shown  here  since  it  is  simllat  to  those  previously 
shown,  the  dynamics  and  feel  system  charact er  ist  ics , as  simulated,  were 
completely  linear.  The  feel  system  dynamics  were  those  of  a conventional 
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(From  Rotoroiu'o  l?1) 


spring-mass-damper  system.  Unfortunately,  Reference  17  does  not  document 

a..  / ^ ; it  was  estimated  as  follows: 

'•p  e 


azr 


(s)  = 


K(s‘  + s + Wj, ) 


^sp  ^ 


2 


U;., 


2Uo  1 
^ x ”^02 


2 


^ZWZ 


1 

T()? 


For  the  six  configurations,  above,  w7  was  sufficiently  greater  than  w 

sp 

that  any  errors  in  the  approximation  of  w or  will  not  significantly 
affect  the  following  results. 


1 


I 


Type  I PTO 

An  appropriate  pilot  model  for  each  of  the  six  configurations  was 
selected  for  the  0 -*  Fg  loop.  From  these,  estimates  of  the  pitch  loop 
resonance  frequencies  w were  made  and  the  phase  criterion  checked  for  the 
a^  -*•  F__  loop  to  ascertain  PIO  susceptibility.  The  results  are  sliown  in 
Table  R.  Two  values  of  <o  are  shown  for  each  case  corresponding  to  a 
nominal  choice  of  pilot  dynamic  equalization  for  optimized  tracking  and  for 
an  unequalized  pilot  (gain  plus  time  delay) — ^ and  respectively. 
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TABLE  8.  ESTIMATED  PIO  PARAMETERS  FOR  CALS PAN  CONE LCURATLONS ; 
TYPE  I PIO 


FLT 

Pilot 

Equal 1 zat ion 

u,rnon 

♦^“Rnom' 

"o 

PIOR 

PIOR 

dust i f led  7 

583 

small  lag 

4 

7.3 

-110 

-210 

4 

Yes 

595 

small  lead 

4.1 

3.8 

-186 

-176 

4.5 

Yes 

596 

moderate  lag 

4.5 

5.0 

-100 

-120 

5 

NO 

601 

large  lead 

4.2 

2.3 

-226 

-141 

4.5 

Yes 

605 

large  lead 

3.9 

3.8 

-170 

-156 

4 

Yes 

610 

lag- lead 

4 

5.5 

-81 

-214 

4 

Yes 

..  - . 



-----  

1(  is  concluded  from  those  data  t Hat  live  ot  tho  six  PIO  cent igurat  ions 
satistv  l l*o  Typo  I PIO  phaso  criterion  ol  present  theory.  Present  theory 
suggests  that  tho  coni Igurat ion  ol  flight  596  is  not  prono  to  tlio  dovolop- 
mont  ol  PIO  as  tho  result  ot  closed  loop  tracking  ol  pitch  attitude. 

Reference  17  notes  the  apparent  contradiction  In* tween  the  Cooper-Harper 
POR  and  the  PIOR  obtained  tor  t light  SSI.  For  this  configuration  the  ratio 

ui  /I  Is  greater  than  one  and  the  Cooper-Harper  rating  ol  pitch  dynamic 
sp  sp  a 

was  Level  1 (POR  » l).  Even  so,  the  PIOR  rating  obtained  was  indicative  ot 
a definite  tendency  to  develop  PIO.  The  pilot's  inexperience  with  the 
rating  scales  and  with  the  flight  test  procedure  was  cited  as  a probable 
reason  tor  this  discrepancy. 

Accotdlng  to  the  present  PIO  theory  there  is  no  reason  to  believe  that 
the  P10R/P0K  results  lor  flight  '>8 ' are  necessarily  contradictory.  It  the 
pilot  rated  only  his  ability  to  control  pitch  attitude  in  the  simulated  IFR 
tracking  task  described  in  Roterenco  17,  then  a Cooper-Harper  rating  ol  l 
might  be  reasonable  (although  there  I n no  way  to  he  certain!;  il  the  pilot 
rated  PIO  tendencies  strictly  according  to  t ho  PIOR  scale  then,  following 
abrupt  eontiol  inputs  with  tight  attitude  control,  present  theory  indicates 


that  unstable  or  nearly  unstable  oscillations  of  the  az^  -*■  F^  loop  can 
develop.  This  result  emphasizes  how  Inadequate  our  understanding  of  the 
basis  for  pilot  rating  really  is;  It  suggests  that  great  care  must  be  taken 
In  assessing  pilot  rating  data  when  more  than  one  cue  or  mode  of  control  may 
exist . 


No  check  was  made  of  the  amplitude  criterion  | az^/ 0( Ju>p jq) I S*  0.012 
since  errors  in  estimation  of  the  acceleration  transfer  function  could 
invalidate  the  result.  Future  investigators  are  encouraged  to  publish  all 
their  data  in  future  experiments  of  these  sort. 

Type  II  PIO 

The  susceptibility  of  the  six  PIO  configurations  to  the  open  loop 
initiation  of  PIO  was  tested  using  the  simplified  criterion  based  on  dominant 
mode  damping  ratio;  viz.,  if  s 0.2,  then  the  configuration  is  suf- 
ficiently resonant  at  frequency  to  initiate  az^  tracking.  If,  further, 
<f>(jw.p)  - -180°  then  the  configuration  is  susceptible  to  Type  II  PIO.  The 
results  are  summarized  in  Table  9.  Note  that  for  these  experiments  the 
stick-free  and  stick-fixed  dynamics  were  the  same  for  the  frequency  range  of 
Interest  here. 

It  is  apparent  that  the  poor  PIOH  obtained  for  flight  59b  could  have 
resulted  from  the  pilot  attempting  to  track  oscillations  of  az  at  frequency 
^ following  control  or  atmospheric  inputs  sufficient  to  excite  the  short- 
When  closed  at  this  frequency,  the  az^  -*•  F^  loop  is 
seen  that  flight  610  is  also  prone  to  Type  II  PIO. 


period  dynamics, 
unstable.  It  is 
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TABLE  l>.  TYPE  II  l’  10  SUSCEPTIBILITY  EVALUATION  OK  CAl.SPAN  CONFIGURATIONS 


Fit 

wsp 

f-sp 

r <9? 

'KJwgp) 

PI  OK 

Potential  lor  Type  11  Pit) 

58J 

8.98 

.66 

No 

n . a . 

4 

Not  possible 

595 

3.72 

.24 

No 

-158 

4.5 

Not  likely 

59b 

6.08 

.13 

Yes 

-200 

5 

Poss ible 

601 

2.  16 

.44 

No 

n.  a . 

4.5 

Not  possible 

605 

3.72 

.19 

Yes 

-148 

4 

Not  likely 

610 

5.47 

.05 

Yes 

-253 

4 

Pons ible 

S umma r v 

The  data  provided  by  Reference  17  appear  to  support  the  present  P10 

theory.  These  same  data  do  not  appear  to  validate  7t  w /I,  (or  its  more 

sp  sp  it 

recent  variations)  as  a reliable  metric  Tor  the  prediction  or  assessment  of 
PtO.  To  this  author,  it  is  clear  that  this  parameter  is  an  indication  ot 
pitch  handling  qualities  but,  according  to  theory  ot  this  report,  this  is 
merely  a necessary  consideration  for  PIO  assessment --not  a sufficient  one. 

What  these  data  do  clearly  Illustrate  Is  how  carefully  we  should  treat 
Cooper-Harper  or  IMOR  scale  data  to  avoid  misconstruing  the  messages  they  may 
provide  due  to  our  lack  of  a unified  understanding  of  the  handling  qualities 
problem  ol  longitudinal  mode  control. 

These  data  appear  to  support  the  simplified  criterion  Cu  - 0.2  as  an 

K 

indication  ot  Type  11  PIO  susceptibility. 
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K. 


A- 7 A (M  = l.l,  h - 15,000  ft) 


* 


The  A- 7 series  of  aircraft  have  no  history  of  longitudinal  1M0  problems. 
It  is  of  interest  to  test  this  airplane  against  the  PIO  theory  to  ascertain 
wltv  PIO  is  not  a problem  (at  least  at  this  one  flight  condition). 

The  (light  condition  selected  for  study  is  not  one  that  is  typical  of 
A-7  operations.  It  was  chosen  because  the  corresponding  airframe  dynamics 
are  similar  to  those  for  the  T-38A.  This  flight  condition  could  be  reached 
as  the  result  of  a supersonic  dive  from  high  altitude. 

With  no  stability  augmentation  the  required  airframe  dynamics  are  (from 
Ko t erence  18): 


0 _ -44. 3(b  + 2.02) 

tSe  ‘S  s[s2  + 2(.  185)  (8.81)s  + 8. HI2] 


rad /rad 


fie  . 222,;[.2  ♦ 2(.072H21,6)s  ft/Mc2/rad 

«e  Is2  + 2 (. 185) (8. 81)s  + 8.812] 

(it  Is  assumed  that  tx  “ 10  ft). 


The  pitch  dynamics  are  shown  on  the  Bode  plot  of  Figure  37;  these  are 
approximately  equal  to  a pure  gain  in  the  region  of  probable  crossover 
frequency.  A reasonable  mode  for  pilot  dynamics  tor  pitch  tracking  would  he 


Yp(s) 


_J^E 

1.0s  + 1 


-.2s 

e 
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The  open-loop,  pilot-airplane  dynamics  are  shown  in  Figure  37  for  the  case 
where  no  control  or  feel  system  dynamics  exist.  The  open-loop  phase  is  also 
shown  for  the  unequalized  pilot.  Without  further  discussion,  it  is  clear 


V 


to 


Figure  37.  A-7A  Pitch  Attitude  Dynamics 

that  the  closed-loop  pitch  dynamics  can  become  resonant  at  *=  4-5.3  rad/ 
sec  with  the  given  pilot  dynamics.  From  the  az  ► dynamics  (Figure  183 
it  is  clear  that  PIO  cannot  develop  over  this  range  of  frequencies.  In  the 
absence  of  control  and  feel  system  dynamics. 

If  the  pilot  decreases  the  level  of  lag  equalization  at  a constant 
static  gain,  then  the  resonant  frequency  can  become  greater  than  5.3  rad/ 
sec.  For  the  unequal i zed  pilot,  the  resonant  I requeney  for  pitch  oscilla- 
tions must  become  greater  than  about  >5  rad/sec  (from  Figure  383  in  order  for 
a PIO  to  develop. 

Resonances  in  pitch  due  to  abrupt  control  or  large  aerodynamic  input 

can  occur  at  the  short-period  frequency  ul  — 8.81.  Since  (,  = 0.185  is 

sp  sp 

v0..’  it  appears  that  these  short -per lod  responses  would  be  subject ivelv 
predictable  (according  to  the  simplified  Type  11  Pit'  criterion).  From 
Figure  18  we  see  that  the  acceleration  loop  dynamics  would  be  unstable  it 
tuned  to  a crossover  I requencv  ot  8.81  rad/sec.  Also, 
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Figure  38.  A-7A  Normal  Acceleration  Dynamics 


7—^(8. 81J  ) 

e 


0.14  >0,012  g/degree/  sec 


Therefore,  conclude  that  the  A-7,  wltli  no  control  or  feel  system  dynamics, 
is  PlO-prone  at  this  flight  condition. 

However,  the  A-7  does  have  significant  control  and  feel  system  dynamics 
even  without  stability  augmentation.  The  airplane  has  two  bobweights  which 
respond  to  0 and  a7  to  create  equivalent  stick  forces;  the  stick-free 
dynamics  are  therefore  substantially  different  from  the  stick-fixed  dynamics 
Figure  Id  illustrates  the  first-approximation  system  schematic.  The  system 
parameters  were  supplied  to  this  author  by  Vought  Aircraft.  When  all  four 
feedback  loops  are  closed,  the  resulting  stick-free  dynamics  at  this  flight 
vendition  become  (using  the  short-hand  notion): 


rad  / lb 


a2p  , . 3501.072.21.6  ] 

Fs  ''S)  ” l . 207 . 7.2811. 448 , 32 . 4 ] 


ft /sec 2 /lb 


The  analysis  of  pilot-vehicle  system  dynamics  will  not  be  repeated;  It  is 
seen  that  t he  bobweight  feedbacks  have  the  principal  effects  of  increasing 
the  short-period  damping  ratio  and  decreasing  the  natural  frequency.  Since 
(.  > 0.2,  and  since  i|>(7.28J)  is  very  nearly  equal  to  -180  degrees,  it  appears 

that  the  probability  for  PIO  initiation  is  significantly  reduced  when 
account  is  taken  of  the  feel  and  control  system  dynamics. 


<r/*‘)  (r  a<l) 


Hobwiltght 


Figure  39.  A-7A  Keel  and  Control  System  Schematic  Description 

The  facts  that  the  stick-free  short-period  damping  ratio  C,  is  0.207, 
and  no  PIO  problems  exist,  are  taken  to  be  further  support  for  the  subjec- 
tive predictability  criterion  - 0.2. 


Olio  final  point  is  worth  noting  about  the  bobweight  effects  on  the  A- 7 
dynamics.  Since  control  system  friction  is  unavoidable,  then  for  small 
amplitude  motions  the  pilot  actually  flys  the  stick-fixed  dynamics.  If  the 
motion  become  large  enough  to  excite  the  bobweights,  then  the  airplane 
dynamics  are  those  of  the  stick-free  responses.  This  change  in  airplane 


dynamics  can  suddenly  occur.  Since  the  stick-free,  short-period  natural 
frequency  is  less  than  the  stick-fixed  value,  then  the  sudden  change  in 
aircraft  dynamics  is  actually  stabilizing  in  the  sense  that,  with  constant 
pilot  gain,  the  closed-loop  system  becomes  more  stable  when  the  bobweight 
feedbacks  are  initiated  following  small  amplitude  motion  tracking.  This 
effect  is  exactly  opposite  that  previously  seen  in  the  T-38A  example  with 
the  original  control  system.  It  is  proposed  that  this  mechanism  provides  a 
reasonable  theoretical  basis  to  support  the  empirically  derived  criterion  of 
Neal  (Reference  22)  tor  the  optimized  design  of  bobweights. 
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SECTION  Vll 


THE  IDENTIFICATION  OF  P 10  IN 
FLIGHT  TEST  OK  SIMULATION 

A.  GROUND-BASED  SIMULATION 

According  to  the  theory  of  this  report , longitudinal,  short-period  PIO 
cannot  exist  in  a fixed-base  simulation;  the  normal  acceleration  cue  is 
essential  to  the  mechanics  oi  PIO.  The  sorts  of  pilot-vehicle  system  oscilla 
tions  that  occur  in  fixed-base  simulations  of  pitch  tracking,  for  example, 
and  which  have  been  called  PIO  in  past  work,  tiave  boon  purposely  excluded 
from  the  PIO  classification  by  the  definition  proposed  in  this  report.  The 
reasons  for  doing  this  have  been  previouslv  discussed. 

However,  the  use  of  fixed-base  simulation  and  the  PIO  theory  may  enable 
the  early  diagnosis  of  PIO.  This  would  require  that  pitch  attitude  tracking 
experiments  be  devised  in  a manner  to  excite  the  extremes  of  piloting 
behavior;  if  it  can  be  determined  by  measurement  that  closed-loop  resonances 
can  occur,  then  an  analytical  check  of  the  acceleration  loop  phase  criterion 

4>(.W  ) - - 1 80° 

K 

at  ouch  resonant  frequency  may  be  sufficient  to  ascertain  whether  or  not  PIO 
would  be  obtained  in  flight.  These  results  from  a combination  of  fixed-base 
simulation  and  from  analysis  could  also  be  used  for  the  preparation  of 
specifications  for  the  motion  drive  system  of  a moving-base  simulation. 

This  might  permit  the  use  of  a moving-base  simulator  to  verify  the  PIO 
tendencies  of  an  aircraft  configuration  during  its  development  and  prior  to 
its  flight. 

A moving-base  simulator  can  only  approximate  the  dynamics  of  flight. 

I'he  chief  limitation  is  the  allowable  range  of  cockpit  travel.  The  usual 
method  for  limiting  the  travel  is  to  introduce  acceleration  washout  in  t he 
mot  ion  drive  system  software  interface  between  the  computed  aircraft  accelera 
t Ion  and  the  servos  and  actuators  that  physical Iv  move  the  cockpit.  Thus, 
the  acceleration  dynamics  become: 
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[r*  (s)]  ■ (8)]  x w(s) 

8 J simulated  8 airplane 


W(s)  • washout  dynamics 


L 


The  typical  washout  system  will  therefore  add  a phase  lead  to  the  pilot-felt 
acceleration  dynamics  that  would  not  actually  exist  In  flight.  According  to 
the  theory  of  this  report,  this  could  eliminate  the  possibility  for  finding 
i'10  in  a moving  base  simulator  (assuming  that  it  existed  in  actual  tlight> 
if,  with  no  washout 

♦(Hio>  v -180° 

and  if,  with  washout 


$(.Uopio)  + 


W(.1«pi0)  > -180° 


Consider  the  YF-17  example  previously  discussed.  P10  was  not  diagnosed 
as  a potential  problem  with  that  aircraft  based  upon  moving-base  simulation 
ot  the  approach  and  landing  problem;  yet  the  in-fright  simulation  ot  the 
airplane  resulted  in  serious  I’lO.  It  has  been  suggested  elsewhere  that  lack 
ot  resolution  of  the  visual  system  presentation  in  the  near-flare  condition 
might  have  invalidated  Northrop's  landing  simulation  of  the  YF-17.  This  is 
conceivable.  It  is  also  conceivable  that  the  "pucker  factor"  was  (rot  suf- 
ficient to  make  the  simulator  pilot  excite  P10.  It  is  the  opinion  ot  this 
author,  however,  that  the  real  problem  was  llkelv  to  have  been  with  the 
motion  drive  svstem. 

Sinacori  (Reference  2 \)  describes  the  mot  ion  drive  logic  in  use  at 
Northrop  at  the  time  ot  t hr-  YF-17  simulation;  he  also  lists  the  mot  ion  drive 
svstem  parameters  appropriate  for  the  YF-17.  It  it  is  assumed  that  these 
are  representative  ot  those  actually  used  tot  the  YF-17  approach  and  landing 
simulation,  then  t ho  normal  acceleration  washout  dynamics  were,  to  a first 
approx imat ion , 
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Till'  effect  of  these  dynamics  on  tlu-  total  phase  ntv’le  of  the  a7  *•  K 

P s 

loop  is  illustrated  in  Figure  40.  The  total  system  nhase  angle  $ is  defined 
as  in  the  previous  YF-17  analysis  (Figure  10):  it  is  assumed  that  t = 0.25. 
Figure  40  clearly  indicates  two  points  of  note: 

o The  acceleration  washout  makes  the  simulated  YF-17  acceleration 
control  loop  significantly  more  stable  titan  is  the  (estimated) 
actual  YF-17  in  the  region  of  probable  ”10  frequencies.  For  w- 
3,  in  fact,  tlie  simulated  YF-17  with  the  original  control  system 
is  not  qualitatively  different  from  the  YF-17  with  the  modified 
control  system. 

o Allowing  for  uncertainties  in  estimates  of  YF-17  dynamics  and  for 
the  actual  motion  drive  dynamics  employed,  it  appears  that  P10  may 
not  have  been  possible  to  produce  in  the  moving  base  simulator — 
regardless  of  "pucker  factor"  or  of  the  fidelity  of  the  visual 
svstem--since  ip  ( 3 j ) is  slightly  greater  than  -180°  witli  washout. 

An  additional  fidelity  problem  could  have  been  created  due  to  the  use 
of  amplitude  sealing  in  the  motion  drive  system: 
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Mi  is  response  rat  io  is  loss  than  tho  or  l tor  ton  value  (.0121  proposed  tor 
sub  loot ivo  prod io t ab 1 1 i t v ot  normal  aooolorat ion  rosponsos. 


In  viow  ot  the  practical  constraints  on  simulation,  it  is  problemat  ioal 
whether  a P 10  problem  can  bo  discovered  in  a simulator  by  an  empirical, 
trial  and  error  approach.  The  use  ot  simulation,  exclusively,  as  a means 
tor  assessing  the  PIO  tendencies  ot  an  aireraf t-oont rol  system  design  is  a 
vet  \ t { s k v approach  to  aircratt  design  and  development. 
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B.  IN-FLIGHT  SIMULATION 


The  use  of  a variable  stability  airplane  for  the  discovery  or  assess- 
ment of  P 10  may  be  an  attractive  alternative  to  moving  base  simulation.  It 
is  not  a panacea,  however. 

It  is,  in  general,  impossible  to  simulate  the  P10  characteristics  of  a 
specif ic  airplane  without  the  use  of  direct  lift  control  (DLC)  in  a variable 
stability  aircraft.  The  use  of  DLC  permits  the  natural  dynamics  of  the 
variable  stability  airplane  to  be  decoupled  in  pitch  and  heave  and,  there- 
fore, to  be  independently  controlled. 


This  was  a potential  error  in  the  simulated  YF-17  flight  tests  with  the 
NT-33A  conducted  prior  to  the  prototype's  first  flight  (Reference  8).  The 
NT-33A  does  not  have  a DLC  system.  In  order  to  match  the  0/6  dynamics  of 


the  YF-17  in  approach  and  landing  it  was  necessary  to  select  the  approach 

airspeed  to  give  the  correct  L(  (or  1/Tg  ) ; this  could  only  be  done  to  a 

first  approximation.  The  parameters  w and  ^ were  matched  using  response 

sp  sp 

feedback.  Since  the  normal  acceleration  dynamics  are  unknown  for  the  actual 
YF-17  and  can  only  be  estimated  for  the  YF-17  as  simulated  by  the  NT-33A,  it 
is  not  possible  to  completely  evaluate  the  degree  to  which  errors  in  simula- 
tion of  az  responses  influenced  the  assessment  of  the  YF-17  handling 
qualities.  It  is  instructive  to  consider  what  these  might  have  been;  con- 
sider the  following: 


1.  For  the  NT-33A  in  the  approach  and  landing  configuration  selected 

to  match  the  YF-17  dynamics, 

U = 140  kt  = 236  ft/sec 
o 

1/T.  =0.9 

0 2 

it  can  be  determined  from  Reference  25  that  for  this  condition: 


=7.39  ft 

X 

Mr  = -4.49 

1/  sec 

Z,  = -15.34 

ft /sec- 
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M - ~ *482  1/sec 
M • » - . 216  1/sue 

ut 

Ma  “ -2. HA  1/sec- 
Za  - -229.4  f c/sec n 

Assume  that  response  feedback  techniques,  only,  were  used  In  the 
YF-17  simulation  of  Reference  8 (i.e.,  model-following,  or  similar 
techniques,  are  assumed  not  to  have  been  used).  It  follows  that 
the  numerators  Ng^(s) , N‘^  (s)  and  N'^zp  (s)  for  the  simulated  YF-17 

dynamics  were  the  same  as  those  for  the  basic  NT-33A . It  Is  noted 
that  this  Is,  In  general,  not  true  when  a ’)1.C  system  is  used. 


For  a rigid  airframe  aZj>  « az-(l.  0.  It  follows  by  direct  synthesis 
that  for  the  YF-17,  as  simulated  hv  the  11T-33A: 


l/p 

17.841.072,7.441 

6 (S) 
e 

1.8 9, 1.98J 

For  the  actual  YF-17  In  the  .approach  and  landing  configuration, 

i - 15.7  ft 
x 

U 118  kt  - 199  ft/sec 
o 

l/Tn  - 0.68 

Up 

The  control  effectiveness  derivatives  are  unknown.  Assume  those  for 
the  F-‘>  in  a similar  flight  condition  are  reasonable  approximations 
(errors  in  these  will  only  affect  the  acceleration  gain). 


Mx  = -5.31  1/sec 

e 

Z\  - 26 .9  t t /sec ' 
e 


*>.  Assume  that  for  tho  Yl‘-17 
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It  has  boon  this  author's  experience 

Mi 

ot  ton  .loouralo  when  V is  nut  '>10. 

* Ge- 
ls an  excel  lout  approximat  ion  to  the 

transfer  Inaction  (I).  It  slum  Id  lu* 

control  system,  in  landing  .’ppro.ich  u 

then  equation  is  invalUI. 


that  this  approximation  is 
Nolo,  also,  that  equal  ion  (2) 
numerator  ot  the  NT- 1 iA 
iiotcil  that  it  the  YF-17  1 J Jght 
sics  I laps  as  an  act  ive  control 


u.  Using  the  approxlm.it  ion  loi  tin'  acco  ler.il  ion  transter  Inaction 
numerator  the  e.-.t  (mated,  a<  . o.i  1 VI  17  transter  I unci  ion  is: 


* *'jp  lti.ti  [ . 080,  i . li  t J 

1 .89,  l . 98 1 


1 1) 


Unfortunately,  there  is  no  wav  to  verity  this  result  without 
more  data. 


It  is  apparent  that  ll  equation  (2)  is  a valid  apt.i ox tmat  ion  lor  the  YF-17 
then  a comparison  ot  equations  (l)  and  (i)  indicates  that  the  NT- l JA  simula- 
tion ot'  the  YK-17  was  possibly  in  error  on  two  counts.  First,  the  magnitude 
ot  static  response  was  about  So  percent  lower  in  the  simulation  than  actual. 
Second,  the  error  in  simulation  produced  more  phase  lag  in  the  simulated 
a7  response  transfer  than  would  actually  exist  in  the  YK-17.  The  numerators 
of  transfer  functions  (1)  and  (3)  are  compared  on  the  Hode  plot  of  Figure  41. 
This  comparison  suggests  that  in  the  vicinity  of  probable  1*  10  frequency  the 
NT-T3A  simulation  of  the  YF-17  may  have  been  slightly  more  l’ 10- prone  than 
the  actual  YF-17.  However,  the  analysis  of  Section  VI  also  concludes  that 
the  YF-17  with  the  original  control  system  was  PlO-prone — therefore  sub- 
stantiating lialspan's  major  conclusion. 
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liK'ite  41.  1‘has**  Anglo  Comparison  of  N *'  ( ) u>)  for  the 

°e 

K- 1 7 In-K  light  Simulation 


C.  Kt.lCHT  I'EST  , ] 

I’.'O  prob!  'ms  with  now  aircralt  designs  otton  1 irst  occur  In  the  latter 
(ages  o t t light  test  or  at  tor  the  aircraft  has  boon  delivered  to  the  opera- 
t ional  congnuuilv.  this  mav  bo  due  to  the  use  ot  t light  tost  techniques 

which  emphasize  Ml  1.-K-878M1  compliance  and  the  determination  ot  the  tradi-  , 

t Ional  measures  ot  airplane  stability  and  control  (o.g.,  short-period  fre- 
quency and  damping  ratio,  stick  force  per  g,  stick  deflection  lor  trim, 

* 

etc.).  When  tests  are  performed  which  emphasize  tracking  U'.g.,  weapons  . 

delivery,  air-to-air  gunnery,  refueling,  etc.)  then  I’lO  may  appear  it  it  is,  i 

indeed,  a latent  problem.  It  would  he  preferable  to  conduct  baseline 

tracking  tests  (e.g.,  as  proposed  in  Reference  .’4)  during  the  initial 

stages  ot  t l ight  testing  t o enable  an  early  diagnosis  of  PIO.  Revised 

flight  test  procedures  and  the  theory  ot  this  report  should  permit  the  rapid 

detection  and  elimination  ot  I’lO  problems. 
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A very  real  problem  that  can  arise  in  flight  testing  is  to  determine 
whether  a control  problem  should  be  called  a PIO  or  a pitch  sensitivity 
problem.  This  question  might  appear  to  be  in  the  same  league  with  deter- 
mining how  many  angels  can  stand  on  the  head  of  a pin.  It  is  not.  If  the 
balance  between  the  pitch  and  heave  degrees  of  freedom  is  understood,  then 
the  opportunities  are  increased  for  eliminating  the  control  problem.  In  an 
era  of  multimode  flight  control  systems,  CCV,  and  direct  force  control,  this 
is  a very  real  benefit  to  the  process  of  flight  control  system  design. 

Once  an  aircraft  reaches  the  flight  test  stage  of  development,  PIO  may 
be  difficult  and  expensive  to  eliminate.  The  normal  acceleration  dynamics 
az^/5i  are  dominated  by  La  and  are  therefore  beyond  control  at  that  point. 

The  elimination  of  PIO,  once  it  is  known  to  exist,  should  probably  be  done 
by 
i 

1.  Reducing  control  system  phase  lags. 

2.  Minimizing  abrupt  changes  in  stick-free  dynamics  due  to  control 
nonlinearities  or  SAS  saturation. 

3.  Improving  pitch  attitude  handling  qualities. 

4.  In  extreme  cases,  adding  a new  dynamic  element  to  the  flight 
control  system  to  attenuate  pitch  attitude  responses. 

Control  stick  dampers  have  been  used  without  (apparently)  great  success  as 
past  "cures"  for  PIO  tendencies;  Lite  documentation  of  these  trials  is  poor, 
it  has  been  suggested  in  past  work  that  the  use  ol  dampers  is  to  be  dis- 
couraged. On  the  basis  of  tills  report,  this  author  is  not  greatly  enthusiastic 
about  a stick  damper,  either.  However,  there  mav  be  cases  where  a damper 
can  be  successfully  used.  An  example  might  be  when  a PIO  results  due  to 

very  small  4 . A damper  cent  igurat  ion  which  decouples  the  feel  system 

sp 

dynamics  might  be  used  to  attenuate  (lie  short-period  response.  This  would 
be  at  tlie  expense  of  Increased  stick  forces. 


UN 


The  use  of  linear  prefilters  to  heavily  attenuate  the  pilot's  (electri- 
cal) output  should  be  discouraged  as  a means  of  preventing  the  transmission 
of  large,  inadvertent  commands  to  surface  actuators  in  a fly-by-wire  system. 
The  control  system  phase  lags  may  lead  to  serious  control  problems  which  can 
escape  diagnosis  in  fixed-  and  moving-base  simulation.  A better  technique 
might  be  to  use  „ mlinear  filtering  logic  which  only  attenuates  the 
amplitudes  of  large,  high-frequency  signals  without  increasing  the  control 
system  phase  lags  within  the  bandwidth  of  the  closed  loop  pitch  control 
loop. 
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SECTION  VIII 

DISCUSSION  AND  CONCLUSIONS 


A PIO  theory  has  been  presented  that  is  believed  to  be  consistent  with 
most  opinions  found  in  the  literature  or  expressed  privately  about  the 
problem.  The  theory  appears  to  unify  divergent  viewpoints  regarding  the 
significance  to  PIO  of  airframe  dynamics,  feel  and  control  dynamics,  task, 
system  nonlinearities,  and  motion  vs  visual  cues.  A number  of  PIO  case 
histories  are  presented  which  appear  to  confirm  the  theory. 

It  is  postulated  that  if  in  the  control  of  pitch  attitude,  the  pilot- 
vehicle  system  dynamics  are  sufficiently  resonant,  and  if  the  pilot-felt 
normal  acceleration  closed-loop  would  be  unstable  when  the  loop  gain  is 
adjusted  to  make  the  crossover  frequency  equal  to  the  resonant  frequency  of 
the  pitch  loop,  then  the  airframe-control  system  dynamics  are  PlO-prone.  A 
PIO  initiated  as  the  result  of  pitch  attitude  tracking  is  defined  as  a 
Type  I PIO. 

Aircraft  configurations  that  feature  a very  lightly  damped  stick-fixed 
or  stick-free  dominant  mode  are  postulated  to  be  susceptible  to  the  initia- 
tion of  PIO  following  abrupt  control  or  atmospheric  disturbance  inputs  of 
magnitude  sufficient  to  excite  the  dominant  mode.  A PIO  initiated  as  the 
result  of  open-loop  control  or  external  inputs  is  defined  as  a Type  II  PIO. 
It  is  proposed  that  a configuration  is  susceptible  to  Type  II  PIO  when  the 
dominant  mode  damping  ratio  is  less  than  0.2. 

The  major  difficulty  in  applying  this  theory  is  to  determine  whether 
the  pitch  loop  dynamics  can  become  sufficiently  resonant  as  the  result  of 
precision,  piloted  control  in  tracking  modes.  The  state  of  the  art  of 
pilot-vehicle  system  dynamics  cannot  entirely  resolve  this  matter.  It  is, 
however,  possible  to  select  models  for  pilot  dynamics  and  perform  a closed 
loop  analysis  to  bound  the  range  of  frequencies  over  which  closed  loop 
resonances  could  occur.  If  the  acceleration  loop  is  unstable  when  closed  at 
any  of  these  frequencies,  then  PIO  is  a possibility.  This,  however,  may  be 
overly  conservative  for  design  purposes.  Validity  probably  must  be 
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determined  on  a case-by-case  basis.  Any  physical  dynamic  system  can  be 
driven  to  resonance  given  sufficient  gain,  for  example;  what  must  be 
assessed  in  a design  or  test  is  the  probability  that  such  resonances  will 
occur  in  actual  operational  usage  of  the  airplane.  This  is  an  engineering 
problem  that  is  probably  best  addressed  from  an  integrated  viewpoint 
blending  ground-based  and  in-flight  simulation,  witli  analysis  used  as  a 
forecasting  tool. 

It  is  noteworthy,  however,  that  subtle  errors  in  the  design  of  flight 
control  systems  have  not  been  responsible  for  those  P 10  cases  that  have  been 
documented.  These  have  resulted  from  deficiencies  that  could,  within  the  . 

present  state  ot  the  art,  now  be  eliminated  witli  t lie  application  of  el  omen- 
tary  pilot-vehicle  system  dynamic  theory.  A principal  source  of  Pit)  has 
been  the  influence  of  artificial  feel  system  dynamics  and  feel  nonlinearities 
(especially  friction)  on  the  dynamic  response  of  the  airplane  to  pilot- 
applied  stick  force.  It  appears  that  pitch  handling  qualities  will  be  poor, 
almost  regardless  of  stick-fixed  dynamics,  when  feel  system  nonlinearities 
make  the  stick-free  dynamics  sensitive  to  amplitude  of  the  motion  such  that 
the  maximum  pilot  gain  required  for  closed  loop  stability  is  less  at  large 
than  at  small  amplitudes.  Problems  of  this  sort  have  most  often  been  due  to 
the  use  of  bobweights  in  past  designs.  The  saturation  of  a rate-  or 

i 

position-limited  SAS  can  produce  the  same  effect. 

Based  upon  the  success  of  this  theory  it  appears  tlvat  the  pilot's 
principal  output  is  stick  force  rather  than  deflection.  This  suggests  two 
corol  lar ies : 

r 

1.  Stick  force  (static  and  dynamic)  should  not  become  "too  small." 

2.  Feel  and  control  system  dynamics  should  not  result  in  stick  * 

deflection  leading  stick  force.  If  that  were  to  happen  then  the 

I 

pilot’s  Internal  mechanisms  for  sensing  and  controlling  force 
output  would  become  unbalanced. 

It  is  conceivable  that  with  small  stick  force  levels  (item  1)  or  with  lead- 
producing  feel  system  dynamics  (item  2)  additional  pilot  lag  is  produced  as 
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a result  of  neuromuscular  system  dynamics.  If  so,  the  effect  would  be  to 
promote  instabilities  or  closed  loop  resonance  in  the  precision  control  of 
pitch  attitude.  A P10  can  then  result  when  the  acceleration  loop  is 
unstable  at  the  resonant  frequency. 

The  importance  of  feel  and  control  system  dynamics  to  the  P10  problem 
is  central  and  cannot  be  overemphasized.  The  present  theory,  however,  can 
provide  no  support  for  the  notion  that  stick  force  per  g is  a PIO  parameter, 
per  se.  As  a general  conclusion,  PIO  tendencies  can  be  minimized  by  mini- 
mizing flight  control  system  phase  lags  in  the  principal  region  of  manual 
control  interest  (nominally  1 < w < 10  rad/sec). 

A major  conclusion  of  this  report  is  that  the  potential  for  PIO  can  be 
predicted  entirely  on  the  basis  of  linear  systems  analysis.  The  basic 
philosophy  is  easily  stated:  there  can  be  no  unwanted,  large  amplitude 
oscillations  so  long  as  the  small  amplitude  motions  are  stable.  The  pre- 
diction of  fully-developed  PIO  frequency  and  amplitude  does,  however,  require 
full  consideration  of  all  major  system  nonlinearities;  even  here,  the  use  of 
simple  quasi-linear  sinusoidal  describing  function  analysis  appears  to  be 
very  successful  as  a diagnostic  tool. 


Appendix 

NORMAL  ACCELERATION  DYNAMICS 


aZp  = aZ  “ ^x^ 


az  = w - UQq 


N*>>  = sN6e(s)  - U°N6e(s) 


N6e(s)  = Z|5e  f®2  " (Mq+Md)s  " (“a  “ z/  Za) 


a a z6eMw_zwM6e 

(s)  = s^xN5  (s)  = s£x(M(Se+Z<seM^)  s + M{  +Z{ 


N/P(s)  = 


aZn 

N5eP^) 


[z6e  - ^x^fie+Z'SeMw)]  s2  + [_Z6e(Mq+M6t)  " *’x(z6eMw-zwM6e)j  s 

/ M<5e  \ 

-Zfie  (* - ^ z*j 

r 1 o z6e(Mq+M(i)  + MZfigMw-ZwMfie) 

[Z<5e  “ 5-x(M6e+Z<5eMv)J  5,2  _ Z{e  - £x (Mfie+Zfie^ 

( M6e  'N 

Z<5e  (^a  “ Z6e  Zoj 
z6e  - ^x(^6e+Z6e^w^ 


This  expression  is  complete  and  assumes  only  that  the  short-period 
approximations  are  valid.  In  generic  form 

NjZp(s)  = Kz(s2  + 2?zwZs  + 


However , 


fHICSWWO  F tat  Ml  AKy 
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Then 


-UoMfi 


e T0„ 


It  often  happens  that  ^xMfig/ZSe  = 2;  assume  this  to  be  true.  Assume  also 
chat  |M«5e|  » |ZgeM^|;  then 


Since,  to  a good  approximation, 

2'-sp“!sp  = " zw  - (Mq+M.) 
it  follows  that 


Z6e 

9 r 4_  ^ I 

+ A1 

KZ 

_-2 Csp  sp  + T0r  ' 

k ZSe  7-1 

or 


It  is  clear  that  is  always  very  small.  Typically,  is  a factor  of  about 
2-5  times  greater  than  uiSp.  The  phase  differences  between  az(t)  and  aZp(t) 
will  therefore  be  significant  at  frequencies  greater  than  u)sp.  At  flight 
conditions  where  ujgp  is  less  than,  or  approximately  equal  to,  the  crossover 
frequency  of  the  pilot-aircraft  system,  the  distinction  between  az  and  aZp 
as  piloting  cues  becomes  of  potential  significance.  Note  that  so  long  as 
r,£  • 0,  aZp(t)  always  leads  az(t)  due  to  the  8 component  of  az,p. 
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The  response  ratio  0/aZp  is  of  probable  importance  to  longitudinal 
handling  qualities.  Using  the  above  approximations,  this  ratio  can  be 


written  as 


Thus,  for  fixed  Cx  and  at  a constant  speed  the  ratio  of  attitude  rate  (or 
attitude)  to  the  pilot's  normal  acceleration  is  entirely  parameterized 
bv  I/Tq  , . 

In  cases  where  i M /Z  =2  is  not  a good  approximation, 

X 0 V’ 


z,  (m  - ~e  z) 

lSe  \ Z>$e  aJ 


Z 7 


'6  - e (M.  +ZC  M-) 

e x 6 6 w' 

e e 


should  be  used. 
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